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Abstract 

Dermal  penetration  of  chemicals  and  drugs  is  important  to  both  toxicologists 
and  pharmacologists.  Drug  developers  try  to  enhance  and  environmental  profession¬ 
als  try  to  limit  penetration  of  chemicals  through  the  skin.  Both  can  use  predictive 
biologically-based  mathematical  models  to  assist  in  understanding  the  processes  in¬ 
volved.  When  these  models  are  based  on  physiological  and  biochemical  parameters 
which  can  be  measured  in  the  laboratory,  they  can  be  extremely  useful.  Appropri¬ 
ately  validated  models  based  on  first  principles  can  be  predictive  of  human  exposures 
when  the  processes  involved  are  adequately  understood.  In  this  thesis  we  develop 
four  new  physiologicaUy-based  pharmacokinetic  (PBPK)  models  to  predict  blood 
concentrations  of  dibromomethane  (DBM)  in  rats  after  neat  hquid  and  vapor  ex¬ 
posure.  These  four  new  models  expand  previously  developed  homogeneous  models 
by  adding  skin  subcompartments.  These  new  models  improve  the  prediction  of  the 
blood  concentrations  especially  early  in  the  exposure.  Sensitivity  analysis  shows  that 
one  of  the  permeabihty  constants  followed  by  the  blood  air  partition  coefficient  have 
the  most  impact  on  blood  concentration  predictions.  With  proper  validation  the 
new  models  could  be  used  to  improve  species,  dose,  and  duration  extrapolations  of 
chemical  or  drug  penetration.  They  could  also  be  used  to  investigate  and  predict 
concentrations  of  drugs  or  chemicals  in  different  parts  of  the  skin. 
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Physiologically-Based  Pharmacokinetic  Modeling 
of  Skin  Absorption  Using  Dermal  Subcompartments 


I.  Introduction 

1.1  Overview 

The  effect  a  chemical  will  have  on  the  human  body  is  usually  unknown  until 
actual  experiments  are  conducted  on  humans.  This  can  be  a  dangerous  proposition 
due  to  the  toxic  hazards  of  some  chemicals,  thus  a  method  of  predicting  these  ef¬ 
fects  would  greatly  decrease  this  danger.  One  method  to  make  these  predictions  is 
to  conduct  experiments  using  laboratory  animals  and  try  to  extrapolate  the  results 
across  species  to  predict  the  effect  on  humans.  Many  researchers  have  developed 
mathematical  models  to  simulate  the  physical  and  chemical  processes  in  the  body. 
A  physiologically-based  pharmacokinetic  (PBPK)  model  is  one  such  model  which 
uses  mathematical  equations  to  represent  the  flow  of  a  chemical  through  the  body. 
Some  applications  of  PBPK  models  include  modeling  inhalation  exposure,  vapor  der¬ 
mal  exposure,  and  liquid  dermal  exposure.  The  equations  of  a  PBPK  model  contain 
physiological  and  biochemical  parameters  that  can  be  measured  in  the  laboratory  or 
estimated.  Using  parameters  that  can  be  measured  in  humans,  as  well  as,  laboratory 
animals  allows  for  the  possibility  of  the  results  to  be  extrapolated  across  species.  It 
can  be  difficult  to  model  dermal  absorption  because  of  the  complex  nature  of  the 
skin  and  its  penetrabihty.  Models  which  include  more  detailed  equations  represent¬ 
ing  the  skin  and  its  subcompartments  can  lead  to  improved  model  predictions  and 
improve  the  possibility  of  extrapolating  the  results.  Sensitivity  analysis  can  identify 
the  parameters  of  a  model  that  have  the  most  significant  impact  on  the  model  out¬ 
put,  which  reveals  the  most  important  parameters  to  measure  in  the  laboratory  or 
estimate  if  necessary. 
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1.2  Problem 


It  is  important  for  the  Air  Force  to  be  able  to  predict  the  effect  of  chemicals 
on  humans,  which  includes  determining  the  concentration  of  the  chemical  in  various 
parts  of  the  body  after  exposure.  Models  have  been  developed  that  predict  blood 
concentrations  in  rats  relatively  well,  however  more  detailed  models  that  use  dermal 
subcompartments  and  are  based  on  biological  parameters  can  give  more  accurate 
predictions  and  make  extrapolation  across  species  to  humans  easier. 

1.3  Scope 

This  research  expands  current  PBPK  models  to  include  dermal  subcompart- 
ments  in  series  and  parallel.  The  models  developed  simulate  percutaneous  absorption 
of  dibromomethane  (DBM)  vapor  and  neat  DBM  liquid.  The  model  determines  the 
concentration  of  the  chemicals  in  various  parts  of  the  body  and  does  not  account 
for  the  physical  change  the  chemical  may  cause  to  the  skin.  Some  of  the  new  model 
parameters  that  were  unable  to  be  measured  had  to  be  optimized.  When  techniques 
to  measure  these  parameter  values  in  the  laboratory  become  available,  the  estimates 
can  be  replaced  with  the  actual  values.  The  accuracy  of  the  models  is  based  on 
comparisons  made  between  the  models^  output  and  experimental  data  collected  in 
the  laboratory  from  experiments  involving  rats.  Metabolism  in  the  Hver  is  the  only 
form  of  metabohsm  that  is  represented  in  this  model.  The  model  results  have  not 
been  extrapolated  across  species  to  predict  effects  on  humans,  however  the  model 
parameters  were  chosen  so  that  future  extrapolation  should  be  possible.  Sensitivity 
analysis  was  performed  to  identify  the  most  important  parameters  to  measure  or 
estimate. 

1.4  Approach 

This  research  created  new  models  that  predict  blood  concentration  profiles  of 
chemicals  after  percutaneous  absorption.  In  this  thesis,  we 
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1.  Reconstructed  a  previously  published  vapor  model  and  verified  the  results  using 
previously  collected  laboratory  animal  data  for  DBM. 

2.  Reconstructed  a  previously  developed  neat  hquid  model  and  verified  the  results 
using  previously  collected  animal  data  for  DBM. 

3.  Reconstructed  a  previously  developed  layered  subcompartment  vapor  model 
and  added  a  few  additional  constraints. 

4.  Extended  the  vapor  model  to  include  a  parallel  follicle  compartment. 

5.  Combined  the  layered  subcompartment  approach  and  the  parallel  subcompart¬ 
ment  approach  into  one  new  model  for  the  vapor  exposure. 

6.  Extended  the  neat  Hquid  model  to  include  separate  sub  compartments  repre¬ 
senting  distinct  layers  of  the  skin. 

7.  Extended  the  neat  Hquid  model  to  include  a  paraUel  foUicle  compartment. 

8.  Combined  the  layered  subcompartment  approach  and  the  paraUel  subcompart¬ 
ment  approach  into  one  new  model  for  the  neat  Hquid  exposure. 

9.  Compared  aU  of  the  models  to  previously  collected  experimental  data  with 
graphical  plots  and  a  calculated  sum  of  squared  deviations. 

10.  Performed  sensitivity  analysis  on  the  models  and  identified  the  most  important 
model  parameters  to  be  measured  or  estimated. 

1.5  Design  Considerations 

The  new  models  developed  in  this  research  should  be  able  to  predict  blood 
concentration  profiles  of  chemicals  other  than  the  ones  discussed  in  this  thesis  by 
using  the  corresponding  biochemical  parameters.  We  chose  DBM  to  use  in  the  vaH- 
dation  of  the  models  because  it  is  one  of  the  chemicals  that  the  Toxicology  Division 
of  the  Armstrong  Laboratory  is  currently  interested  in  modeHng  [24].  The  new  pa¬ 
rameters  derived  in  the  new  models  are  biological-based  so  the  parameters  should 
be  measurable  in  humans  so  the  results  can  be  extrapolated  across  species. 
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1.6  Summary  of  Thesis 

The  thesis  is  organized  as  follows: 

Chapter  11  reviews  the  fundamentals  of  skin  absorption  and  physiologically- 
based  pharmacokinetic  (PBPK)  modeling.  It  also  includes  a  discussion  of  the  previ¬ 
ous  and  current  research  being  done  in  this  area. 

Chapter  III  presents  the  model  design  and  implementation. 

Chapter  IV  presents  the  results  of  the  new  models  compared  to  animal  data 
previously  collected  in  the  laboratory.  It  also  includes  sensitivity  analysis  which 
identifies  the  parameters  that  have  the  largest  impact  on  the  model  output. 

Chapter  V  presents  a  summary  of  the  work  completed  and  gives  the  conclusions 
reached.  It  also  includes  recommendations  for  future  research. 
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IL  Background 


2.1  Overview 

Both  pharmacologists  and  toxicologists  are  concerned  with  the  absorption  of 
chemicals  through  the  skin.  The  ability  to  predict  the  concentration  of  a  chemical 
in  various  parts  of  the  body  following  percutaneous  absorption,  which  is  the  absorp¬ 
tion  of  chemicals  through  the  skin,  can  be  obtained  with  a  good  understanding  the 
absorption  process.  To  fully  understand  percutaneous  absorption  one  must  first  un¬ 
derstand  the  complex  structure  and  function  of  the  skin.  There  are  many  things  that 
affect  the  permeabiUty  of  the  skin  including  regional  and  structural  variations  and 
chemical  properties.  Pharmacokinetic  models  can  be  developed  to  simulate  the  phar¬ 
macokinetics  (absorption,  distribution,  metabolism,  and  elimination)  of  a  chemical 
in  the  body. 

Not  all  models  use  the  same  approach  to  represent  diffusion  through  the  skin. 
Models  that  include  biologically-based  parameters  that  can  be  measured  in  different 
species  have  the  best  chance  for  extrapolation.  Models  based  on  mathematics  alone 
can  be  very  accurate  for  the  specific  conditions  for  which  they  are  created,  but  extrap¬ 
olation  to  other  conditions  and  species  can  be  extremely  difficult.  Physiologically- 
based  pharmacokinetic  (PBPK)  models  are  based  on  physiological  and  biological 
parameters  that  can  be  measured  in  different  species.  Sensitivity  analysis  of  PBPK 
models  can  be  used  to  identify  the  most  important  biologically-based  parameters  to 
measure. 

Before  explaining  our  model  design  and  implementation,  we  wiU  review  the 
structure  of  mammalian  skin,  some  previous  work,  PBPK  modeling,  a  previously 
developed  dermal  vapor  absorption  model,  and  a  previously  developed  dermal  layered 
subcompartment  model. 
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2.2  Skin  Structure 


The  skin  is  a  very  important  part  of  our  body.  The  skin  is  the  largest  organ 
in  the  body  by  weight  and  is  the  primary  means  of  identifying  individuals.  The 
entire  body  is  enveloped  by  the  skin  allowing  it  to  protect  the  body  from  the  outside 
environment  which  is  its  most  important  function.  Bodily  fluids  would  escape  and 
external  fluids  would  enter  the  body  much  easier  if  the  skin  was  not  present.  An 
animal’s  skin  enables  it  in  certain  ways  to  communicate,  respond,  and  adapt  to  its 
environment.  The  skin  can  protect,  camouflage,  warn,  or  attract  attention,  and 
in  some  animals  its  glands  can  secrete  substances  that  can  repel  or  attract  other 
animals  [26]. 

The  skin  has  many  detailed  structural  and  functional  differences  among  various 
species  but  the  overall  effect  of  the  skin  on  the  body  is  the  same  in  all  species  [26]. 
Certain  characteristics  of  the  skin  vary  over  different  regions  of  the  body  implying 
the  skin  is  not  a  uniform  organ.  For  example,  the  thickness  and  permeability  of  the 
skin  differs  greatly  from  some  areas  of  the  body  to  others  [46] .  The  skin  is  as  thick 
as  400-600  [i  in  plantar  and  palmer  callus  versus  10-20  /jl  thick  on  the  back,  arms, 
legs,  and  abdomen  [30].  The  number  of  skin  appendages  and  other  structures  also 
vary  across  the  body. 

2.2.1  Skin  Layers.  Most  people  normally  think  of  the  skin  as  a  single¬ 
layered  membrane  that  covers  and  protects  the  body.  A  more  detailed  description 
identifies  two  distinct  but  hnked  layers  which  contain  sublayers  within  them.  The 
outer  layer  is  the  epidermis  and  the  inner  layer  is  known  as  the  dermis  or  cornium 
[33].  Figure  2.1  shows  a  cross  section  of  the  skin  revealing  its  complex  structure. 

The  epidermis  contains  two  parts,  the  stratum  corneum  and  the  viable  epider¬ 
mis.  The  stratum  corneum  is  the  outermost  part  of  the  epidermis  and  often  referred 
to  as  the  horny  layer  of  the  skin.  The  stratum  corneum  is  made  up  of  layers  of 
non-living  cells  which  are  constantly  being  replaced.  These  non-living  cells  originate 
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Figure  2.1  Cross  section  view  of  the  skin 
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in  the  germinative  basal  layer  of  the  epidermis  and  gradually  move  upwards  replen¬ 
ishing  the  stratum  corneum  about  every  two  to  four  weeks  [25,  33].  The  stratum 
corneum  is  the  main  barrier  to  penetration  through  the  skin  [36].  The  viable  epider¬ 
mis  not  only  suppHes  the  dead  cells  to  the  stratum  corneum  but  it  also  surrounds 
hair  follicles,  even  into  lower  skin  layers.  The  epidermal  layer  makes  up  about  25 
percent  of  the  total  skin  thickness  [32]. 

The  dermis  or  cornium  is  the  inner  layer  of  the  skin  which  accounts  for  the 
bulk  of  the  skin  thickness  and  about  15  to  20  percent  of  the  total  body  weight.  The 
dermis  is  a  matrix  of  loose  connective  tissue  which  is  traversed  by  blood  vessels, 
nerves,  and  lymphatics.  The  dermis  is  also  penetrated  from  above  by  hair  follicles 
and  sweat  glands.  This  tough,  resilient  and  viscoelastic  tissue  has  many  functions 
which  include  providing  nourishment  to  the  epidermis  [26]. 

2.2.2  Skin  Appendages.  A  relatively  small  fractional  area,  approximately 
10“^  in  humans,  of  the  skin  is  covered  with  skin  appendages.  Skin  appendages 
include  hair  follicles,  sweat  ducts,  and  sebaceous  glands  [36].  These  appendages  help 
to  regulate  body  temperature  by  sweating  through  the  sweat  glands.  The  hair  that 
grows  from  hair  follicles  helps  keep  the  body  warm  by  providing  a  layer  of  insulation. 
The  number  of  skin  appendages  varies  from  one  species  to  another  and  from  one 
region  of  the  body  to  another.  The  skin  appendages  extend  from  the  surface  of  the 
skin  into  and  through  the  dermis  where  a  dense  capillary  mesh  surrounds  them  [35] . 

2.2.3  Skin  Permeability.  Skin  permeabiUty  is  important  to  understand 
when  concerned  with  skin  absorption  because  it  determines  how  much  the  body 
accumulates  or  eliminates  substances  through  the  skin.  Skin  permeability  varies 
across  the  body  with  the  regional  variations  in  the  physiological  properties  of  the 
skin  [38].  The  skin  permeabihty  is  also  different  in  each  layer  of  the  skin,  the  stratum 
corneum  being  the  least  permeable  layer.  In  places  where  the  skin  is  diseased, 
damaged,  or  removed,  absorption  of  water-soluble  chemicals  can  be  1,000  times 
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greater  than  normal.  Even  through  normal  intact  skin,  there  can  be  a  huge  difference 
in  the  absorption  of  chemicals  with  different  properties  [36].  These  variances  in  skin 
permeabihty  make  modeling  skin  absorption  very  challenging. 


2.3  Classical  Pharmacokinetic  Modeling 

Laboratory  animals  are  often  used  in  experiments  to  determine  the  effect  of 
percutaneously  absorbed  chemicals.  Chemical  concentration  in  blood  taken  from 
the  animals  at  set  times  after  exposure  gives  some  indication  of  the  processes  taking 
place.  These  in  vivo  techniques  give  a  more  complete  analysis  than  in  vitro  tech¬ 
niques  since  metaboUsm,  nervous,  and  humoral  responses  are  not  present  in  in  vitro 
experiments  [21].  It  is  important  to  be  able  to  develop  models  that  make  predictions 
consistent  with  the  data  from  in  vivo  experiments  so  fewer  laboratory  animals  will 
have  to  be  sacrificed.  These  models  can  be  used  to  predict  the  effect  of  chemicals  on 
humans. 

There  are  two  basic  types  of  models  that  have  been  used  for  modehng  percuta¬ 
neous  absorption,  one  is  the  compartment  model  [12,  45]  which  is  based  on  first-order 
rate  constants  and  the  other  is  the  diffusion  model  which  is  based  on  the  diffusion 
laws.  The  diffusion  model  has  been  used  most  often  when  modehng  absorption  of 
solutes  through  the  skin  and  it  is  the  most  correct  mechanistically  [40]. 

The  mathematics  used  to  describe  skin  permeabihty  in  the  diffusion  model  are 
based  on  Pick’s  second  law  which  is  presented  in  Equation  (2.1). 


dC  _  d^C 
dt  ~^dx^ 

Where 

C  =  Concentration 
t  =  Time 


(2.1) 
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D  =  Average  membrane  diffusion  coefficient  (distance^  time 
X  =  Distance 

The  following  integrated  form  of  Pick’s  law  represents  the  steady-state  flux  of 
solute  [36]. 


,  D-{Ct-C2) 

8 

Where 

Jg  =  Steady-state  flux  of  solute  (moles  distance"^  tirne"^) 

8  =  membrane  thickness  (cm) 

It  is  also  necessary  to  include  a  term  to  describe  the  relationship  between 
external  and  surface  concentrations  which  can  be  done  using  the  solvent-membrane 
distribution  coefficient  {Km)-  Including  this  coefficient,  Equation  (2.2)  becomes  the 
following. 


,  Km-D-ACg 
8 

Where 

ACa  —  Concentration  difference  of  solute  across  membrane  (moles  cm“^) 

Km  =  (Solute  sorbed  per  cc  of  tissue)/(Solute  in  solution  per  cc  of  solvent) 

The  permeability  coefficient  (fep)  is  described  using  a  portion  of  this  flux  equa¬ 
tion  and  is  normally  used  in  a  model.  The  permeabihty  coefficient  describes  the  rate 
of  absorption  of  the  chemical  (distance/time). 


(2.3) 


kp 


Km-D 


8 


Where 


(2.4) 
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kp  =  Permeability  constant  for  solute  (distance  time 

The  effective  permeability  of  each  skin  layer  relative  to  the  permeability  of 
the  whole  skin  can  be  obtained.  The  following  equation  shows  the  relationship  be¬ 
tween  the  permeability  constants  for  the  individual  skin  layers  and  the  permeability 
constant  for  the  whole  skin  [5,  34], 


1 


Where 


(2.5) 


”  PermeabiKty  constant  of  the  whole  skin 
kp.  =  PermeabiKty  constant  for  the  layer  of  the  skin 

Several  mathematical  models  have  been  developed  that  represent  the  absorp¬ 
tion  of  chemicals  through  the  skin  using  a  layered  subcompartment  approach.  Some 
of  the  layered  models  use  two  [11,  27,  41,  45]  or  three  [40]  subcompartments. 

The  skin  appendages  act  as  diffusional  shunts  for  penetration  through  the  skin. 
Solutions  penetrating  the  skin  through  the  skin  appendages  by-pass  the  upper  layers 
of  the  skin,  thus  exhibiting  different  permeabiKty  characteristics  from  the  rest  of  the 
skin.  The  effect  of  parallel  pathways  on  permeabiKty  relative  to  the  whole  skin  can 
theoreticaUy  be  obtained.  The  foUowing  equation  shows  the  relationship  between 
permeabiKty  constants  for  the  parallel  pathways  and  the  permeabiKty  constant  for 
the  whole  skin  [7]. 


K.,  ■  =  E(*b..  ■  Ai)  (2-6) 

i 

Where 

Aak  =  Area  of  the  whole  skin 

kp-  =  Permeability  constant  for  the  parallel  pathway  of  the  skin 
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Ai  =  Area  of  the  parallel  pathway  of  the  skin 

Many  researchers  have  also  investigated  the  effect  of  skin  appendages  as  a  shunt 
pathway  in  percutaneous  absorption  [14,  17,  35,  15,  43,  44].  Although  researchers 
agree  that  appendages  act  as  a  diffusional  shunt  through  the  skin,  the  impact  of  ap¬ 
pendages  on  percutaneous  absorption  varies,  and  is  related  to  the  intrinsic  physical 
properties  of  the  chemical  tested  and  the  time  of  apphcation  [14].  During  the  initial 
stage  of  absorption,  the  appendages  may  have  a  significant  contribution,  however 
once  the  chemical  penetrates  the  stratum  corneum  and  the  viable  epidermis,  the 
contribution  of  the  appendages  to  the  overall  absorption  process  is  insignificant  [17]. 
Experiments  have  revealed  that  differences  in  skin  structure,  for  example  skin  ap¬ 
pendages,  could  have  an  effect  on  percutaneous  absorption  when  a  slowly  absorbed 
test  penetrant  is  used  [37].  Some  mathematical  models  have  been  developed  which 
contain  a  parallel  skin  subcompartment  to  represent  a  shunt  pathway  for  chemicals 
penetrating  the  skin  through  skin  appendages  [2,  4,  18,  45]. 

Different  approaches  have  been  taken  to  represent  the  diffusion  equations 
within  the  context  of  a  model.  Some  models  use  partial  differential  equations,  solv¬ 
ing  them  with  Laplace  transformations  [1,  10,  27,  47].  This  approach  allows  the  skin 
concentration  to  vary  with  skin  depth  and  with  time.  These  partial  differential  equa¬ 
tions  make  it  difficult  to  use  biologically-based  parameters  and  they  have  boundary 
conditions  that  can  be  difficult  to  incorporate  into  a  model. 

Most  of  the  previously  developed  models  mentioned  above  are  primarily  de¬ 
scriptive  because  they  are  based  on  classical  pharmacokinetic  techniques  and  other 
techniques  that  are  not  based  on  biological  parameters  which  can  be  measured  or 
estimated  for  the  purpose  of  extrapolation  from  one  species  to  another.  Most  of  the 
time,  the  equations  in  these  models  are  fit  to  a  curve  of  the  experimental  data  and 
are  only  useful  for  certain  chemicals,  species,  and  exposure  concentrations.  Since 
these  models  are  based  on  mathematical  analysis  alone,  they  do  not  give  a  good 
biological  description  of  the  skin  or  the  processes  that  take  place  within  it.  Another 
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type  of  model  has  been  developed  that  is  based  more  on  physiological  properties  but 
contains  the  mathematics  to  describe  the  diffusion  process  accurately. 

2.J^  Physiologically-Based  Pharmacokinetic  Modeling 

Using  the  experimental  data  to  predict  the  effect  of  the  chemical  on  other 
species  is  challenging.  One  tool  used  to  make  these  predictions  is  a  mathematical 
model  which  includes  physiological  and  pharmacokinetic  principles.  A  physiologically- 
based  pharmacokinetic  (PBPK)  model  is  one  such  model  which  uses  mathematical 
equations  to  describe  the  dynamics  of  chemicals  moving  through  the  body.  A  PBPK 
model  can  represent  some  of  the  nonlinear  processes  in  the  body  by  using  physi¬ 
ological  parameters  such  as  blood  flow  rates,  breathing  rates,  blood  volumes,  tis¬ 
sue  volumes,  permeability  of  membranes,  and  partitioning  of  chemicals  into  tissues 
[21,  42]. 

Some  of  the  parameters  in  a  PBPK  model  can  be  measured  in  a  laboratory, 
however  some  of  the  parameters  must  be  estimated  to  complete  the  model.  One 
technique  for  estimating  the  unknown  parameters  is  to  optimize  the  unknown  pa¬ 
rameters  in  the  model  to  get  the  best  fit  to  the  corresponding  experimental  data. 
These  estimated  values  can  be  used  until  actual  measurements  can  be  made  in  the 
laboratory. 

2.4.1  Mass-Balance  Equations.  PBPK  models  represent  the  body  by  iden¬ 
tifying  tissue  compartments  made  up  of  single  organs  or  groups  of  physiologically 
similar  organs  and  tissues.  This  body  representation  is  an  improvement  over  the 
classical  pharmacokinetic  models  which  are  based  on  first  order  transfer  rates  and 
represent  the  body  by  identifying  one  to  three  compartments  based  on  a  semilogarith- 
mic  plot  of  plasma  concentration  versus  time  [9].  A  PBPK  model  uses  mass-balance 
equations  to  represent  the  flux  through  a  tissue  compartment  and  the  processes  that 
take  place  within  it.  Mass  balance  equations  can  be  used  to  ensure  that  the  total 
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mass  of  the  chemical  in  the  body  is  being  conserved  throughout  the  tissue  compart¬ 
ments.  The  following  equation  is  a  general  mass-balance  equation  representing  a 
tissue  compartment. 


=  Q- 
dt 


(2.7) 


In  Equation  (2.7),  i  refers  to  tissue  compartment  i,  V  refers  to  volume,  Q  refers 
to  the  blood  flow,  C  refers  to  the  concentration,  and  Ri/i,  refers  to  the  partition 
coefficient  between  the  tissue  and  blood  [21,  29].  A  partition  coefficient  is  the  ratio 
of  the  distribution  of  a  chemical  between  body  components  at  steady-state  [23]. 
Some  of  the  tissue  compartments  will  have  more  complex  mass-balance  equations 
to  represent  other  chemical  activity  within  the  compartment.  The  other  chemical 
activities  represented  could  be  air-to-blood  exchange  in  the  lungs,  metabolism  in  the 
liver,  or  absorption  from  the  surface  into  the  skin.  When  the  system  of  mass-balance 
equations  is  solved  simultaneously,  the  concentration  of  the  chemical  in  each  tissue 
compartment  can  be  obtained,  as  well  as,  the  concentration  in  the  blood. 


2. 4-2  Modeling  Using  Dermal  Subcompartments.  PBPK  models  that  iden¬ 
tify  dermal  subcompartments  can  more  accurately  represent  the  physiology  of  the 
skin  and  thus  improve  model  predictions.  Subcompartment  models  can  include 
biologically-based  parameters  which  can  be  measured  in  different  species  and  will 
allow  extrapolation  across  species  to  be  easier. 

Using  subcompartments  in  series  to  represent  the  individual  layers  of  the  skin 
can  be  more  accurate  than  using  a  homogeneous  skin  compartment  since  the  layers 
have  very  different  permeabihty  properties.  A  parallel  subcompartment  to  represent 
the  shunt  pathway  of  skin  appendages  can  also  create  a  more  accurate  model  since 
the  appendages  have  been  shown  to  have  an  effect  on  skin  absorption. 
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When  incorporating  dermal  sub  compartments  into  a  PBPK  model  it  is  neces¬ 
sary  to  created  new  mass-balance  equations  to  represent  the  skin  compartment.  The 
new  mass-balance  equations  introduce  new  parameters  such  as  permeabihties,  parti¬ 
tion  coefficients  between  the  subcompartments,  fractional  areas,  fractional  volumes, 
and  fractional  blood  flows.  Some  of  these  new  parameters  can  be  determined  experi¬ 
mentally  in  the  laboratory,  however  the  unknown  parameters  must  be  estimated  until 
measurements  techniques  are  available.  When  estimating  the  unknown  parameters, 
care  must  be  taken  to  preserve  certain  relationships,  for  instance  the  permeability 
relationships  in  Equations  (2.5)  and  (2.6).  A  relationship  between  the  partition  co¬ 
efficients  of  the  sub  compartments  and  the  partition  coefiicient  of  the  skin  must  be 
maintained  when  modeling  dermal  subcompartments  presented  in  Equation  (2.8). 

Rsklair  •  '  K)  (2-8) 

i 

2,4 ’3  Extrapolation  Across  Species.  The  overall  goal  of  a  PBPK  model 
is  to  be  able  to  extrapolate  the  results  of  the  model  across  species  to  predict  the 
effect  of  a  chemical  on  humans  [21].  The  properties  and  methods  used  in  inter¬ 
species  extrapolation  include  scaling  physiological  parameters,  physiological  time, 
and  physiologically-based  pharmacokinetics.  This  extrapolation  is  not  easy  due  to 
the  physiological  differences  between  laboratory  animals  and  humans.  Therefore,  a 
PBPK  model  based  on  physiological  parameters  that  can  be  measured  in  humans, 
as  well  as,  laboratory  animals  can  make  the  extrapolation  process  much  easier  [31]. 
In  the  absence  of  species  specific  parameters,  body  weight  scahng  or  surface  area 
scaUng  can  be  used.  These  scahng  techniques  depend  on  the  characteristics  of  the 
chemical  being  used  [42]. 

2.4 -4  Sensitivity  Analysis.  Measurement  and  estimation  of  biochemical 
and  metabohc  parameters  is  one  reason  extrapolation  of  model  results  across  species 
is  difficult.  Sensitivity  analysis  can  help  determine  which  of  the  parameters  are  the 
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most  important  to  the  model’s  output.  This  will  reveal  which  parameters  should  be 
measured  or  estimated  most  carefully  [13]. 

There  are  several  ways  to  perform  sensitivity  analysis  on  model  parameters 
and  the  best  method  for  a  model  depends  on  the  information  that  is  needed  about 
the  parameters  and  the  model.  One  way  to  determine  the  sensitivity  of  a  parameter 
is  to  observe  the  percent  change  in  the  model’s  output  resulting  from  a  certain 
percent  change  in  the  parameter  [39].  One  method  of  sensitivity  analysis  that  uses 
this  approach  involves  calculating  a  ”log-normahzed  sensitivity  coefficient”.  This 
coefficient  is  a  ratio  of  the  resulting  percent  change  in  the  model  output  to  the 
percent  change  in  the  parameter  [6].  A  Monte  Carlo  simulation  using  a  Latin- 
hypercube  procedure  which  produces  estimates  of  the  differential  sensitivities  of  a 
model  is  sometimes  used  in  sensitivity  analysis  [8] . 

The  sensitivity  of  a  parameter  will  sometimes  change  over  time  and  with  dif¬ 
ferent  exposure  concentrations.  Some  parameters  may  have  an  inverse  affect  on  the 
output  of  a  model.  That  is,  an  increase  in  the  parameter  value  will  cause  a  decrease 
in  the  model’s  output.  AU  of  these  considerations  should  be  looked  at  before  making 
any  conclusions  about  the  sensitivity  of  parameters. 

2.5  A  Previously  Developed  PBPK  Model  for  Dermal  Vapor  Exposure 

McDougal  et  al.  (1986)  developed  a  PBPK  model  for  dermal  absorption  of 
vapors  in  the  rat.  Their  model  was  developed  to  predict  the  blood  concentrations 
of  dihalomethanes  after  dermal  exposure  to  chemical  vapors.  Figure  2.2  shows  the 
lumped  anatomical  compartments  of  their  model.  This  body  representation  was 
patterned  after  a  model  developed  by  Ramsey  and  Anderson  (1984)  which  predicted 
the  blood  concentration  after  inhalation  of  styrene  vapor  in  rats  .  These  compart¬ 
ments  were  originally  chosen  because  they  allow  the  representation  of  metabolism 
and  other  processes  that  axe  important  to  the  distribution  of  chemicals  in  the  body. 
In  this  model,  the  skin  is  represented  by  a  single  homogeneous  compartment  (Fig- 
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Figure  2.3  McDougal  et  aVs  Representation  of  the  Skin  Compartment 


Slowly  Perfused  Tissue: 


v^^  =  Q 


dt 
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sp/b  j 


(2.11) 


Liver: 


kfCl-Vl  Vmaa:  '  Q 
Ri/b  K  •  Ri/b  +  Cl 


(2.12) 


Fat: 


-5= a 


(2.13) 


The  following  are  the  mass-balance  equations  their  model  used  for  the  arterial 
blood  and  the  venous  blood  returning  to  the  lungs. 


Arterial: 


Ga  = 


Qc  •  Cv 


(2.14) 
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Venous: 


(2.15) 


C,= 


Qc 


Where: 
a  —  arterial 
A  =  area  (distance^) 
al  =  alveolar 
b  =  blood 

C  =  concentration  (mass/volume) 
c  =  cardiac 
/  =  fat 

i  =  tissue  compartment  i 

K  =  Michaelis  constant  -  metabolism  (mass/volume) 
kf  =  first-order  metaboHc  rate  constant  (time~^) 

I  =  liver 

P  =  permeability  (distance /time) 
p  =  pulmonary 
Q  =  flow  (volume/time) 

R  =  partition  coefficient  (ratio  of  concentrations) 
rp  =  rapidly  perfused 
sfc  =  surface 
sk  =  skin 

sp  =  slowly  perfused 
t  =  time 
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V  =  volume 


V  =  venous 

Vmax  =  maodmum  reaction  velocity  (mass/time) 

The  parameters  for  their  model  were  obtained  from  the  hterature  or  measured 
experimentally  in  the  laboratory  and  are  presented  in  Section  4.2.1.  Their  model 
was  vahdated  for  several  different  dihalomethanes  and  at  several  different  exposure 
concentrations. 

As  opposed  to  using  a  partial  differential  equation,  McDougal  et  al.’s  model 
uses  an  ordinary  differential  equation  to  represent  the  diffusion  equation  within  the 
context  of  the  model.  The  ordinary  differential  equation  approach  assumes  the  skin 
compartment  is  well-stirred  and  thus  the  change  in  chemical  concentration  within  the 
skin  changes  with  time  and  not  with  distance  through  the  skin.  This  representation 
of  the  diffusion  equation  allows  the  model  to  incorporate  parameters  that  are  based 
on  biological  properties  (blood  flow,  tissue  volumes,  body  weight,  etc.)  of  the  animal 
and  physiochemical  properties  (permeabihty  constants,  partition  coefficients,  etc.) 
of  the  chemical  being  used  in  the  model,  as  opposed  to  using  multiple  first-order 
rate  constants.  The  advantage  of  using  these  types  of  parameters  is  that  they  can 
be  measured  (or  optimized  if  absolutely  necessary)  in  the  laboratory.  This  allows 
the  ability  to  extrapolate  the  results  of  the  models  to  other  species  after  appropriate 
parameters  are  measured. 

S.6  A  Previously  Developed  Layered  Dermal  Subcompartment  Model 

McDaniel  (1993)  extended  the  vapor  model  of  McDougal  et  al.  (1986)  by 
adding  multiple  layered  dermal  subcompartments.  His  models  use  the  same  mass- 
balance  equations  as  McDougal  et  al.’s  vapor  model  except  for  the  equation  for  the 
skin  compartment  (Equation  (2.9)).  McDaniel  developed  three  new  models  that  rep¬ 
resented  the  skin  as  a  series  of  subcompartments  and  compared  the  model  predictions 
to  experimental  data. 
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His  first  model  was  a  two  subcompartment  model  that  identified  a  stratum 
corneum  subcompartment  and  combines  the  viable  epidermis,  dermis,  and  subcuta¬ 
neous  fat  into  a  single  composite  dermal  subcompartment  which  is  where  the  blood 
exchange  takes  place  (Figure  2.4).  This  model  assumes  the  volume  of  the  composite 
dermal  subcompartment  is  the  same  as  the  volume  of  the  stratum  corneum  subcom¬ 
partment.  Equations  (2.16)  and  (2.17)  are  the  new  mass-balance  equations  for  the 
skin  sub  compartment.  This  model  made  blood  concentration  predictions  that  were 
more  accurate  than  the  homogeneous  skin  model. 


blood  flow 


Surface 


Stratum  Corneum 


Composite  Dermal 


Figure  2.4  Two  Layered  Subcompartment  Model 


Stratum  Corneum  Sub  compartment; 


dt 


-  P 8C  '  ^8C  I  Gsfc 


a. 


R 


8c/8fc > 


+  Pcd  ■  ^cd  1  Ccd 


(2.16) 


Composite  Dermal  Subcompartment: 


Vcd — 7^  =  Pcd  •  Acd  (  r> - "b  ^cd 


dt 


Rs 


i/cd 


(2.17) 


Where 
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sc  —  stratum  corneum 


cd  =  composite  dermal 

McDaniel’s  second  model  used  the  same  equations  as  his  first  model  to  rep¬ 
resent  two  skin  subcompartments,  however  in  the  second  model  the  volume  of  the 
composite  dermal  subcompartment  is  assumed  to  be  ten  times  greater  than  the  vol¬ 
ume  of  the  stratum  corneum  subcompartment.  This  model  is  slightly  less  accurate 
in  predicting  blood  concentrations  than  the  first  model  but  is  stiU  a  significant  im¬ 
provement  over  the  homogeneous  model,  as  well  as,  being  more  biologically  correct. 

McDaniel’s  third  model  represents  the  skin  as  three  layered  subcompartments. 
The  three  subcompartment  model  (Figure  2.5)  identifies  a  stratum  corneum  subcom¬ 
partment  hke  the  two  subcompartment  model,  however  the  three  sub  compartment 
model  breaks  the  composite  dermal  subcompartment  into  a  viable  epidermis  subcom¬ 
partment  and  a  new  composite  dermal  subcompartment  consisting  of  the  remaining 
layers  of  the  skin  (dermis  and  subcutaneous  fat).  The  new  composite  dermal  subcom¬ 
partment  is  the  only  skin  subcompartment  in  which  blood  exchange  takes  place.  This 
more  detailed  representation  improved  the  representation  of  the  physiology  of  the 
skin.  The  volume  of  the  viable  epidermis  subcompartment  is  assumed  to  be  three 
times  greater  than  the  volume  of  the  stratum  corneum  sub  compartment  and  the 
volume  of  the  new  composite  dermal  subcompartment  is  assumed  to  be  seven  times 
greater  than  the  volume  of  the  stratum  corneum  subcompartment.  This  assumption 
preserves  the  size  of  the  stratum  corneum  of  the  second  two  subcompartment  model. 
The  mass-balance  equations  for  the  three  skin  subcompartments  are  shown  below. 

Stratum  Corneum  Subcompartment: 


dt 


—  P  -  A 

—  ^  SC 


a 


a 


sfc 


R. 


-f-  P„e  •  ve  ae- 


a 


sc/sfc J 


R. 


scjve  ^ 


(2.18) 
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Figure  2.5  Three  Layered  Subcompartment  Model 


Viable  Epidermis  Subcompartment: 


dCve  _  p  . 

yye  ve  ^ve 
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Cve  I  +  Pcd  *  ^cd  (  Ccd 
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■sc/ve 


Pve/cd 
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Composite  Dermal  Subcompartment: 


dt 


^  Rve/cd 


-cjj 


+  Qcd  (  Ca  “ 


Cca 


R. 


•Cdjhj 


(2.19) 


(2.20) 


Where 

ue  =  viable  epidermis 

The  three  sub  compartment  model  is  the  most  accurate  of  the  three  models  for 
predicting  blood  concentrations  and  is  the  most  biologically  descriptive.  The  new 
models  introduced  several  new  parameters,  some  of  which  can  be  measured  (vol¬ 
umes)  and  others  which  had  to  be  optimized  (permeabihty  constants  and  partition 
coefficients).  The  unknown  parameters  were  optimized  by  a  linear  search  technique 
where  the  parameters  were  varied  one  at  a  time  to  determine  the  value  which  led  to 
the  smallest  error  between  the  model  predictions  and  the  experimental  data.  Mc- 
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Daniel’s  models  were  a  significant  improvement  over  the  homogeneous  skin  model  in 
predicting  blood  concentrations  and  in  representing  the  physiology  of  the  skin. 

2.1  A  Previously  Developed  PBPK  Model  for  Dermal  Neat  Liquid  Exposure 

McDougal  et  al.  developed  a  PBPK  model  to  represent  the  absorption  of  a  neat 
Hquid  chemical  after  dermal  exposure.  This  model  is  similar  to  the  vapor  exposure 
model  in  that  it  uses  the  same  body  compartment  representation  (Figure  2.2).  The 
hquid  model  has  one  additional  mass-balance  equation  to  represent  the  amount  of 
the  chemical  on  the  surface  of  the  skin  (Equation  (2.21)). 

Amount  on  Skin  Surface: 

^=P.k-  A.,  ( -  C.,,)  (2.21) 

dt  \-^8k/lq  j 


Where 

Iq  =  Liquid 

In  this  model  the  area  of  the  skin  (A,*)  is  the  area  exposed  to  the  chemical 
and  not  the  entire  surface  area  of  the  body  as  in  the  vapor  model.  The  hquid 
model  also  uses  different  equations  for  the  concentration  in  the  venous  blood  and 
the  concentration  in  the  arterial  blood.  The  vapor  model  assumed  these  blood 
concentrations  reached  steady-state  so  quickly  compared  to  the  concentrations  in 
the  other  body  compartments  that  the  rate  of  change  in  concentration  is  equal 
to  zero,  thus  simphfying  the  mass-balance  equations  (Equations  (2.14)  and  (2.15)). 
Equations  (2.22)  and  (2.23)  are  the  new  equations  for  the  arterial  blood  concentration 
and  the  venous  blood  concentration. 

Arterial: 

=  (2.22) 

dt  y  ^b/air  / 
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Venous: 


(2.23) 


K 


dt 


-  {Qc  ■  C.) 


The  parameters  for  their  model  were  obtained  from  the  literature  or  measured 
experimentally  in  the  laboratory  and  are  presented  in  Section  4.4.1.  Their  model 
has  been  validated  using  several  chemicals  by  comparing  the  model  results  to  ex¬ 
perimental  data  collected  from  rats  in  the  laboratory.  This  model,  like  the  vapor 
model,  contains  biologically-based  parameters  that  can  be  measured  or  estimated 
for  extrapolation  across  species. 


2. 8  Summary 

This  review  uncovered  some  of  the  complex  structures  and  functions  of  the  skin, 
as  well  as,  some  of  the  factors  that  influence  the  absorption  of  chemicals  through  it. 
The  review  of  the  previous  work  revealed  the  basic  approaches  used  in  pharmacoki¬ 
netic  modehng  of  percutaneous  absorption  and  some  of  their  drawbacks.  A  general 
discussion  of  PBPK  modeling  provided  some  of  the  background  necessary  to  under¬ 
stand  how  it  improves  upon  previous  techniques  of  pharmacokinetic  modehng.  We 
used  the  background  discussed  in  this  chapter  as  the  foundation  for  the  new  models 
we  developed  in  this  research. 
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III.  Model  Design  and  Implementation 

3. 1  Overview 

Some  Physiologically  Based  Pharmacokinetic  (PBPK)  models  have  been  de¬ 
veloped  to  predict  the  effect  of  chemicals  in  a  rat  after  dermal  absorption.  These 
models  are  based  on  biological  and  physiological  parameters  that  can  be  measured 
or  estimated  so  the  results  can  be  extrapolated  across  species. 

We  have  taken  these  previously  developed  models  and  extended  them  to  in¬ 
clude  dermal  subcompartments  which  represent  some  of  the  layers  of  the  skin  and 
the  skin  appendages.  The  dermal  subcompartments  give  a  more  detailed  biological 
description  of  the  skin  and  predict  blood  concentrations  more  accurately. 

3.2  Computational  Methods 

To  solve  the  systems  of  ordinary  differential  mass-balance  equations,  we  used  a 
modified  Fortran  77  program  originally  written  by  Dr.  Quinn.  This  program  solved 
the  system  of  equations  by  calling  the  subroutine  IVPAG  from  the  International 
Mathematic  and  Statistics  Library  (IMSL).  The  subroutine  IVPAG  was  created  to 
solve  an  initial  value  problem  for  ordinary  differential  equations  using  an  Adams- 
Moulton  or  Gear  Method.  We  also  converted  some  of  our  models  to  ACSL,  which 
is  the  programing  language  the  Toxicology  Division  of  the  Armstrong  Laboratory 
uses  for  their  skin  absorption  models,  to  ensure  our  method  of  solving  the  system  of 
equations  was  consistent  with  theirs. 

3.3  Optimization  Technique 

Our  program  also  included  a  subroutine  containing  the  method  of  steepest 
descent  and  Newton’s  Method  algorithms  which  were  used  to  optimize  the  unknown 
parameters  to  achieve  the  smallest  weighted  sum  of  squared  deviations  between  the 
model  predictions  and  the  experimental  data  (Equation  (3.1))  [3,  28]. 
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SSD  =  Sum  of  squared  deviations 
n  =  Number  of  experimental  data  points 

The  weight  for  each  exposure  was  taken  to  be  the  average  of  all  of  the  experi¬ 
mental  data  points  for  that  exposure  (Equation  (3.2)).  A  weighted  sum  of  squared 
deviations  was  used  to  ensure  the  optimization  techniques  did  not  favor  a  specific 
exposure  concentration  or  time. 


3-4  Homogeneous  Vapor  Model 

We  reconstructed  McDougal  et  aVs  homogeneous  vapor  model  presented  in 
Section  2.5  and  used  it  as  a  base  to  compare  the  predictions  of  the  new  subcompart¬ 
ment  vapor  models.  The  only  thing  that  changed  as  we  developed  the  subcompart¬ 
ment  models  was  the  equation  for  the  skin  compartment. 


3.5  Layered  Suhcompartment  Vapor  Model 

The  layered  subcompartment  vapor  models  are  based  on  McDaniel’s  models 
presented  in  Section  2.6.  The  layered  sub  compartment  models  used  in  this  research 
have  the  same  set  of  mass-balance  equations  as  McDaniel’s  models.  The  permeability 
constraints  of  Equations  (2.5)  and  (2.6)  were  included.  We  constrained  the  partition 
coefficients  for  the  layered  subcompartments  according  to  Equation  (2.8).  Instead  of 
the  1:10  and  the  1:3:7  estimated  ratios  that  McDaniel  used,  we  used  new  measured 
values  for  the  depths  of  the  layered  subcompartments.  We  used  a  stratum  corneum 
depth  of  11/560  the  depth  of  the  skin,  a  viable  epidermis  depth  of  22/560  the  depth 
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of  the  skin,  and  a  composite  dermal  depth  of  527/560  the  depth  of  the  skin  [20].  Also, 
a  conservation  of  mass  calculations  was  added  to  ensure  that  all  of  the  chemical  was 
accounted  for  at  all  times  throughout  the  simulation  (Equation  (3.3)). 


TABS  —  AMsk  +  AMf  +  AMi  -f  AMgp  +  AM^p  +  AM^x  +  AMmet  (^•^) 
Where 

TABS  =  Total  amount  absorbed  (mass) 

AM  =  Amount  (mass) 
ex  =  Exhaled 
met  =  MetaboUzed 

For  all  of  the  layered  subcompartment  models  used  in  this  research,  the  areas 
of  the  layered  subcompartments  (Age,  .d„e,  and  A^d)  are  assumed  to  be  the  same  as 
the  overall  area  of  the  skin  (A,fc).  Also,  the  blood  flow  to  the  composite  dermal  sub¬ 
compartments  (Qcd)  is  assumed  to  be  the  same  as  the  blood  flow  to  the  overall  skin 
(Qak)  because  the  composite  dermal  subcompartment  is  the  only  subcompartment 
in  which  blood  exchange  takes  place. 

3.5.1  Two  Layered  Subcompartment  Vapor  Model.  The  two  layered  sub¬ 
compartment  model  introduced  four  new  parameters  that  had  to  be  measured  or 
optimized.  The  new  parameters  were  the  stratum  corneum  air  partition  coefficient 
(Rsc/air),  composite  dermal  air  partition  coefficient  (Rcd/air),  the  permeability  con¬ 
stant  for  the  stratum  corneum  subcompartment  (Psc),  and  the  permeability  constant 
for  the  composite  dermal  subcompartment  (Pad).  The  permeability  constant  for  the 
composite  dermal  subcompartment  was  not  optimized  directly  because  it  was  deter¬ 
mined  at  run  time  by  the  permeabihty  constant  for  the  stratum  corneum  subcom¬ 
partment  through  the  following  relationship  which  is  based  on  Equation  (2.5): 
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(3.4) 


1  _  1  J_ 

P,k  Psc  ^  Pcd 

Thus  the  equation  for  Pcd  is: 

The  partition  coefficient  for  the  composite  dermal  subcompartment  was  not 
optimized  directly  because  it  was  determined  at  run  time  by  the  partition  coefficient 
for  the  stratum  corneum  subcompartment  through  the  following  relationship  which 
is  based  on  Equation  (2.8): 


Psk/air  ’  ^sk  —  Psc/air  "  "b  Pcd/air  ’  ^cd 


Thus  the  equation  for  Rcd/air  is: 


(3.6) 


R. 


R. 


cdjair  — 


kja 


V>k  -  Rs 


V.C 


Vcd 


(3.7) 


Because  of  the  uncertainty  as  to  which  subcompartment  permeabihty  and  par¬ 
tition  coefficient  to  optimize  directly  and  which  to  optimize  indirectly  using  Equa¬ 
tions  (3.4)  and  (3.6),  we  also  optimized  Pcd  and  Rcd/air  directly  and  allowed  P,c 
and  Rsc/air  to  be  optimized  indirectly.  This  process  verified  that  the  optimization 
procedure  would  produce  the  same  results  using  both  approaches. 


3.5.2  Three  Layered  Subcompartment  Vapor  Model.  The  three  layered 
sub  compartment  model  introduced  six  new  parameters  that  had  to  be  measured  or 
optimized.  The  new  parameters  were  the  stratum  corneum  air  partition  coefficient 
{Rsc/air),  viable  epidermis  air  partition  coefficient  {Rue /air),  composite  dermal  air 
partition  coefficient  {Rcd/air),  the  permeabihty  constant  for  the  stratum  corneum 
subcompartment  {Psc),  the  permeabihty  constant  for  the  viable  epidermis  subcom- 
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partment  (Pve),  and  the  permeability  constant  for  the  composite  dermal  subcom¬ 
partment  (Pcd)-  The  values  of  i2,c/air  and  P,c  can  be  taken  from  the  two  layered 
subcompartment  model.  We  could  not  use  the  value  of  Rcdfair  Pcd  from  the  two 
layered  subcompartment  model  because  the  composite  dermal  sub  compartment  for 
the  two  layered  subcompartment  model  includes  the  viable  epidermis  and  the  com¬ 
posite  dermal  subcompartment  for  the  three  layered  subcompartment  model  does 
not.  The  permeabihty  constant  for  the  composite  dermal  subcompartment  was  not 
optimized  directly  because  it  was  determined  at  run  time  by  the  permeability  con¬ 
stant  for  the  stratum  corneum  and  viable  epidermis  subcompartments  through  the 
following  relationship  which  is  based  on  Equation  (2.5): 


1 

Psk 


1 

Pcd 


Thus  the  equation  for  Pcd  is: 


(3.8) 


Pcd  = 


Psk  •  Psc  •  Pv< 


Pve  —  Psk  ■  Pve  —  Psk  '  Psi 


(3.9) 


The  partition  coefficient  for  the  composite  dermal  subcompartment  was  not 
optimized  directly  because  it  was  determined  at  run  time  by  the  partition  coeffi¬ 
cients  for  the  stratum  corneum  and  viable  epidermis  subcompartments  through  the 
following  relationship  which  is  based  on  Equation  (2.8): 


Psk/air  '  ^sk  —  Psc/ air  '  ^sc  Pve/air  *  hye  “t”  Pcd/ air  ’  "^cd  (3.10) 


Thus  the  equation  for  Pcd/air  is: 


Pcd/ air  — 


Psk/air  "  ^sk  Psc/ air  ’  ^sc  Pve/air  ‘ 

V^d 


(3.11) 
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As  in  the  two  layered  subcompartment  model,  we  also  optimized  Pcd  and 
Rcd/aiT  directly  and  allowed  and  i?t,e/air  to  be  optimized  indirectly  using  Equa¬ 
tions  (3.8)  and  (3.10).  This  method,  as  opposed  to  the  one  mentioned  above,  ver¬ 
ified  that  the  optimization  procedure  would  produce  the  same  results  using  both 
approaches. 


3. 6  Parallel  Subcompartment  Vapor  Model 

We  developed  a  parallel  subcompartment  vapor  model  to  represent  the  shunt 
pathway  of  skin  appendages  through  the  skin.  This  model  expanded  the  homoge¬ 
neous  vapor  model  by  representing  the  skin  as  two  parallel  subcompartments  (Fig¬ 
ure  3.1).  The  parallel  subcompartment  model  identifies  a  follicle  sub  compartment 
and  a  composite  dermal  subcompartment.  The  foUicle  subcompartment  includes 
the  hair  foUicles  and  sweat  glands  and  the  composite  dermal  subcompartment  in¬ 
cludes  the  normal  layers  of  the  skin  (stratum  corneum,  viable  epidermis,  dermis, 
and  subcutaneous  fat).  The  composite  dermal  subcompartment  for  the  parallel  sub¬ 
compartment  model  is  different  from  the  composite  dermal  subcompartments  of  the 
layered  subcompartment  models.  The  fractioncd  area  of  the  follicle  subcompartment 
{Afo)  is  one  one-hundredth  the  area  of  the  entire  skin  (A,*)  [20].  Blood  exchange 
takes  place  in  both  subcompartments  and  the  fractional  blood  flow  to  the  follicle 
subcompartment  {Qfo)  is  one  fourth  the  blood  flow  to  the  entire  skin  {Qtk)  [20]. 
The  depth  of  the  follicle  subcompartment  is  assumed  to  be  388/560  the  depth  of  the 
entire  skin  [20].  The  following  are  the  mass-balance  equations  for  the  two  parallel 
sub  comp  art  ment  s . 

FoUicle  Subcompartment: 


T/ 


—  Pfo*  ^fo  I  ^sfc 


Cio 


R 


•fo/sfc 


+  Qfi 


("-t 


fo_ 

76, 


(3.12) 
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Figure  3.1  Parallel  Subcompartment  Vapor  Model 


Composite  Dermal  Sub  compartment: 


14,^  =  -P=,  •  f  C.f,  -  .5^1  +Q^(c.-^)  (3.13) 

dt  y  -^cd/ afc  /  \  -^cd/b  / 


Where 
fo  =  Follicle 

Equations  (3.12)  and  (3.13)  replace  Equation  (2.9)  in  the  homogeneous  model. 
The  parallel  subcompartment  model  introduced  four  new  parameters  that  had  to  be 
measured  or  optimized.  The  new  parameters  are  the  follicle  air  partition  coefficient 
{Rfo/air)>  the  composite  dermal  air  partition  coefficient  (Rcd/air)}  the  follicle  perme¬ 
ability  constant  (F/o)>  ^■iid  the  composite  dermal  permeabihty  constant  {Pcd)-  The 
permeability  constant  for  the  composite  dermal  subcompartment  was  not  optimized 
directly  because  it  was  determined  at  run  time  by  the  permeability  constant  for  the 
follicle  subcompartment  through  the  following  relationship  which  is  an  application 
of  Equation  (2.6): 


Psk  ’  -Asfe  =  Pfo  ’  Afo  -|-  Pcd  ‘  Acd 


(3.M) 


3-7 


Thus  the  equation  for  Pcd  is: 


Pcd  = 


P >k  '  ^ak  Pfo  '  -^fo 
Acd 


(3.16) 


The  partition  coefficient  for  the  composite  dermal  subcompartment  was  not 
optimized  directly  because  it  was  determined  at  run  time  by  the  partition  coefficient 
for  the  follicle  subcompartment  through  the  following  relationship  which  is  based  on 
Equation  (2.8): 


Psk/air  '  ^sk  —  Pfo/air  *  "b  Pcd/air  *  "^cd 


(3.16) 


Thus  the  equation  for  Pcd/air  is: 


P. 


P. 


cdlair 


\k/air  *  ^sk  Rfo/air  ’ 


(3.17) 


As  in  the  layered  subcompartment  models,  we  also  optimized  Pcd  and  Pcd/air  di¬ 
rectly  and  allowed  Pfo  and  Pfo/air  to  be  optimized  indirectly  using  Equations  (3.14)  and  (3.16). 
This  method,  as  opposed  to  the  one  mentioned  above,  verified  that  the  optimization 
procedure  would  produce  the  same  results  using  both  approaches. 


3. 7  Parallel-Layered  Subcompartment  Vapor  Model 

We  developed  a  new  model  that  combined  the  two  layered  subcompartment 
model  and  the  parallel  subcompartment  model.  This  new  parallel-layered  subcom¬ 
partment  model  expanded  the  homogeneous  vapor  model  by  representing  the  skin 
as  three  subcompartments  which  are  the  follicle,  stratum  corneum  and  the  com¬ 
posite  dermal  (Figure  3.2).  The  folhcle  subcompartment  includes  the  hair  follicles 
and  sweat  glands.  The  stratum  corneum  subcompartment  consists  of  the  stratum 
corneum  and  the  composite  dermal  subcompartment  includes  the  remaining  layers 
of  the  skin  (viable  epidermis,  dermis,  and  subcutaneous  fat).  The  composite  dermal 
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subcompartment  for  the  parallel-layered  subcompartment  model  is  different  from  the 
previously  discussed  composite  dermal  subcompartments.  The  fractional  area  of  the 
follicle  subcompartment  (Afo)  is  one  one-hundredth  the  area  of  the  entire  skin  {A,k) 
[20].  Blood  exchange  takes  place  in  both  sub  compartments  and  the  fractional  blood 
flow  to  the  follicle  subcompartment  (Qfo)  is  one  fourth  the  blood  flow  to  the  entire 
skin  (Qsk)  [20].  The  depth  of  the  folhcle  subcompartment  is  assumed  to  be  388/560 
the  depth  of  the  entire  skin  [20].  The  depth  of  the  composite  dermal  subcompart¬ 
ment  is  assumed  to  be  549/560  the  depth  of  the  entire  skin.  The  areas  of  the  layered 
subcompartments  {Age  Emd  A^d)  are  assumed  to  be  the  same  as  the  overall  area  of 
the  non-folhcle  portion  of  the  skin  {Acd  from  the  parallel  subcompartment  model). 
Also,  the  blood  flow  to  the  composite  dermal  subcompartments  {Qcd)  is  assumed 
to  be  the  same  as  the  blood  flow  to  the  overall  non-folhcle  portion  of  the  skin  {Qcd 
from  the  parallel  subcompartment  model)  because  the  composite  dermal  subcom¬ 
partment  is  the  only  layered  subcompartment  in  which  blood  exchange  takes  place. 
The  following  are  the  mass-balance  equations  for  the  three  skin  sub  compartments  of 
the  parallel-layered  subcompartment  model. 


Surface  . 

_ _ L 

Follicle 

_ LL_ 

T 

Stratum  Com 

eum 

Composite  Di 
_ 1 

ermal 

J _ 

arterial  blood 


venous  blood 


Figure  3.2  Parallel-Layered  Subcompartment  Model 
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Follicle  Subcompartment: 


dC 


So 


dt 


Pfo  '  ^fo  I  ^sfc 


Cl,  \ 


Pfo/sfcJ 


+  Qfi 


/bl 


(3.18) 


Stratum  Corneum  Subcompartment: 


dCsc  _  p  .  (n_ 

^8C  —  ^8C  '  ^8C  \^afc 


acjsjc  j 


+  Pcd  •  Acd  (  Ccd  — 


a 


R. 


'Bc/cd  y 


(3.19) 


Composite  Dermal  Subcompartment: 
dCcd 


Vcd- 


dt 


Pcd  '  -^cd  ( - C'cd  ]  +  Qcd 


R 


acjcd 


(^“  Rcdibj 


(3.20) 


Equations  (3.18),  (3.19),  and  (3.20)  replace  Equation  (2.9)  in  the  homogeneous 
model.  The  parallel-layered  subcompartment  model  introduced  six  new  parameters 
that  had  to  be  measured  or  optimized.  The  new  parameters  are  the  follicle  air  parti¬ 
tion  coefficient  (ii/o/air),  the  stratum  corneum  air  partition  coefficient  {Rscjair)^  the 
composite  dermal  air  partition  coefficient  {Rcdlair)-}  the  foUicle  permeability  constant 
{Pfo)^  the  stratum  corneum  permeability  constant  (Psc)j  s^^id  the  composite  der¬ 
mal  permeability  constant  {Pcd)’  The  value  of  Rfo/air  3'iid  Pfo  were  taken  from  the 
parallel  subcompartment  model  since  the  follicle  subcompartment  is  physiologically 
the  same  for  both  models.  The  parameters  for  the  stratum  corneum  subcompart¬ 
ment  and  the  composite  dermal  subcompartment  had  to  be  optimized  because  these 
subcompartments  are  different  from  the  ones  in  the  two  layered  subcompartment 
model  because  in  this  model  they  do  not  include  the  follicle  region  of  the  skin.  The 
permeabiUty  constant  for  the  composite  dermal  subcompartment  was  not  optimized 
directly  because  it  was  determined  at  run  time  by  the  permeability  constant  for  the 
stratum  corneum  subcompartment  through  the  following  relationship  which  is  an 
application  of  Equations  (2.5)  and  (2.6). 
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Psk  •  Ask  =  Pfo  •  Afo  + 


1  1 
+ 


-1 


led 


5c  Pcd{Fig  3,2)/ 

Notice  that  Equation  (3.21)  is  the  same  as  Equation  (3.14)  except  that 


(3.21) 


Pcd(Fig  3.1) 


J_  +  1 

Psc  P :d{Fig  3.2) 


-1 


(3.22) 


where  cd{Fig  3.1)  is  the  composite  dermal  sub  compartment  of  Figure  3.1  and 
cd{Fig  3.2)  is  the  composite  dermal  subcompartment  of  Figure  3.2.  Thus  the  equa¬ 
tion  for  Pcd{Fig  3.2)  is: 


p  _  Psc  ’  Psk  ’  Agk  Psc  *  Pfo  ‘  Afo 

“  Psc'Aca-Psk-A^k  +  Pfo^Afo 


(3.23) 


The  partition  coefficient  for  the  composite  dermal  subcompartment  was  not 
optimized  directly  because  it  was  determined  at  run  time  by  the  partition  coeffi¬ 
cients  for  the  stratum  corneum  and  follicle  subcompartments  through  the  following 
relationship  which  is  based  on  Equation  (2.8): 


hjfc  Psc/air  '  hjc  H”  Ped/air  ■  Ved  F  Pfo/ air  ’  h/o 


(3.24) 


Thus  the  equation  for  Ped/air  is: 


P. 


'cdjaiT 


_  Pakjair  *  ^ak  Rac/air  ’  ^ac  Rfo/air  '  ^fo 


Ved 


(3.25) 


We  also  optimized  Ped  and  Ped/air  directly  and  allowed  Pgc  and  Pac/air  to  be 
optimized  indirectly  using  Equations  (3.21)  and  (3.24).  This  method,  as  opposed  to 
the  one  mentioned  above,  verified  that  the  optimization  procedure  would  produce 
the  same  results  using  both  approaches. 


3-11 


3.8  Homogeneous  Liquid  Model 

We  reconstructed  McDougal  et  aVs  homogeneous  liquid  model  presented  in 
Section  2.7  and  used  it  as  a  base  to  compare  the  predictions  of  the  subcompartment 
Hquid  models.  The  mass-balance  equations  for  these  subcompartment  models  are 
the  same  as  the  ones  in  McDougal  et  o/.’s  model  except  for  the  equation  for  the  skin 
compartment. 

3.9  Liquid  Subcompartment  Models 

We  used  the  same  approach  as  we  did  with  the  homogeneous  vapor  model 
to  develop  four  new  hquid  subcompartment  models  (two  layered  subcompartment 
model,  three  layered  subcompartment  model,  parallel  subcompartment  model,  and 
parallel-layered  subcompartment  model).  The  equations  for  the  skin  subcompart¬ 
ments  in  each  hquid  subcompartment  model  are  almost  the  same  as  the  ones  in 
the  corresponding  vapor  subcompartment  models.  The  only  difference  is  that  in 
the  hquid  models  we  use  subcompartment  hquid  partition  coefficients  instead  of 
subcompartment  surface  partition  coefficients. 

3.10  Sensitivity  Analysis 

It  is  important  to  know  which  parameters  have  the  largest  effect  on  the  models’ 
blood  concentration  predictions.  We  performed  sensitivity  analysis  to  identify  the 
skin  compartment  parameters  that  cause  the  largest  percent  change  in  arterial  blood 
concentration  predictions  when  they  are  changed  by  a  set  percent  from  their  original 
values.  One  at  a  time,  we  increased  and  decreased  the  values  of  the  parameters  by 
five  percent  and  computed  blood  concentration  profiles.  We  chose  a  time  early  in 
the  exposure  and  a  time  late  in  the  exposure  and  determined  the  amount  the  output 
changed  due  to  a  five  percent  change  in  the  parameters. 
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3.11  Summary 

We  have  taken  the  previously  developed  homogeneous  models  and  extended 
them  to  include  dermal  subcompartments  in  series  to  represent  the  layers  of  the 
skin  and  in  parallel  to  represent  the  skin  appendages.  Although  the  new  parameters 
are  biologically  based,  some  of  them  had  to  be  estimated  until  better  measurement 
techniques  are  developed  in  the  laboratory.  The  new  sub  compartment  models  pre¬ 
dict  blood  concentrations  more  accurately  because  they  provide  a  better  biological 
description  of  the  skin  and  they  have  more  parameters  that  are  optimized. 
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IV.  Model  Results 


4.1  Overview 

The  predictions  of  the  new  models  were  compared  to  experimental  data  in 
which  rats  were  exposed  certain  chemicals  for  specified  periods  of  time  which  the 
models  are  set  up  to  simulate.  The  sum  of  squared  deviations  and  graphical  plots 
are  used  to  measure  how  well  the  models’  predictions  match  the  experimental  data. 
The  subcompartment  models’  results  were  compared  to  the  homogeneous  models’ 
results  to  ensure  the  new  subcompartment  models  were  at  least  as  accurate  as  the 
original  models.  Sensitivity  analysis  performed  on  the  skin  compartment  parameters 
in  the  models  showed  which  of  the  parameters  have  the  largest  effect  on  the  predicted 
arterial  blood  concentrations  of  the  chemicals. 

4.2  Vapor  Model  Predictions 

The  blood  concentration  predictions  from  the  five  vapor  models  were  compared 
to  previously  collected  data  in  which  rats  were  exposed  to  dibromomethane  (DBM) 
vapor  at  concentrations  of  500  ppm,  1,000  ppm,  5,000  ppm,  and  10,000  ppm  for 
four  hour  periods  [22].  The  increase  in  accuracy  from  the  homogeneous  model  to 
the  subcompartment  models  may  not  be  immediately  obvious  since  the  graphs  are 
displayed  with  a  logarithmic  scale  and  the  relative  differences  between  the  actual 
data  and  the  predictions  are  already  small. 

4.2.1  Homogeneous  Vapor  Model.  The  parameters  for  the  homogeneous 
model  were  determined  experimentally  or  found  in  the  literature  (Tables  4.1  and  4.2) 
[23].  Figure  4.1  shows  how  the  predictions  of  the  previously  published  homogeneous 
model  compare  to  the  experimental  data.  The  homogeneous  model  predicted  the 
concentrations  in  the  blood  reasonably  well  having  a  sum  of  squared  deviations  of 
0.55  when  compared  to  the  experimental  data. 
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Compartment 

Tissue/air  partition 
coefficient 

Blood  flow 
(  %  cardiac  output) 

Volume 

(  %  body  weight) 

Skin 

266.0 

5 

10 

Rapidly  Perfused 

68.1 

56 

5 

Slowly  Perfused 

40.5 

10 

65 

Fat 

792.0 

9 

7 

Liver 

68.1 

20 

4 

Blood 

74.1 

N/A 

N/A 

Table  4.1  DBM  Parameters  for  the  tissue  compartments  of  the  homogenous  vapor 
model 


Constants 

Values 

Vmax  (mg/hr/kg) 

12.5 

Km  (mg/Uter) 

0.4 

Kf  (hr-^-kg-^) 

3.4 

Kp  (cm/hr) 

1.32 

Table  4.2  DBM  Constants  for  the  homogenous  vapor  model 

^.2.2  Two  Layered  Subcompartment  Vapor  Model.  The  unknown  sub- 
compartment  permeabihty  constants  and  partition  coefficients  of  the  two  layered 
subcompartment  model  were  estimated  and  the  optimized  parameter  values  can  be 
seen  in  Table  4.3.  The  two  layered  subcompartment  model  matched  the  data  better 
than  the  homogeneous  model,  yielding  a  sum  of  squared  deviations  of  0.40  when 
compared  to  the  experimental  data.  The  two  layered  subcompartment  model’s  sum 
of  squared  deviations  is  a  27  percent  improvement  in  accuracy  over  the  homogeneous 
model.  Figure  4.2  shows  how  well  the  predicted  values  of  the  two  layered  subcom¬ 
partment  model  match  the  experimental  data.  The  two  layered  sub  compartment 
model  predicts  the  concentration  of  the  chemical  in  the  blood  more  accurately  than 
the  homogeneous  model,  because  it  is  a  more  physiologically  detailed  representation 
of  the  skin  and  it  has  two  more  parameters  that  were  optimized. 
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concentration  (mg/L) 


1000 


DBM  Homogenous  Vapor  Mode! 


Figure  4.1  DBM  Homogeneous  Vapor  Model  Results.  This  graph  shows  how  the 
model’s  predictions  compare  to  the  experimental  data. 


Unknown 

Parameters 

Optimized 

Values 

Psc 

1.46 

Pcd 

14.00 

^8c/ air 

1530.52 

Rcdj  air 

240.66 

Sum  of  Squared 
Deviations 

0.40 

Table  4.3  DBM  Two  Layered  Subcompartment  Vapor  Model  Estimated 
Parameters 
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DBM  Two  Layered  Subcompartment  Vapor  Model 


Figure  4.2  DBM  Two  Layered  Subcompartment  Vapor  Model  Results.  This  graph 
shows  how  the  model’s  predictions  compare  to  the  experimental  data. 

4-2.3  Three  Layered  Subcompartment  Vapor  Model.  The  unknown  sub¬ 
compartment  permeability  constants  and  partition  coefficients  of  the  three  layered 
subcompartment  model  were  estimated  and  the  optimized  parameter  values  can  be 
seen  in  Table  4.4.  The  three  layered  subcompartment  model  yielded  about  the 
same  results  as  the  two  layered  subcompartment  model  yielding  a  sum  of  squared 
deviations  of  0.40  when  compared  to  the  experimental  data.  The  three  layered  sub¬ 
compartment  model’s  sum  of  squared  deviations  is  a  27  percent  improvement  in 
accuracy  over  the  homogeneous  model.  Figure  4.3  shows  how  well  the  predicted 
concentration  values  of  the  three  layered  subcompartment  model  match  the  experi¬ 
mental  data.  The  three  layered  subcompartment  model  yields  more  accurate  blood 
concentration  predictions,  as  well  as,  a  more  detailed  representation  of  the  structure 
of  the  skin. 
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Unknown 

Optimized 

Parameters 

Values 

Pec 

1.46 

Pve 

17.47 

Pcd 

70.36 

Rsc/air 

1530.53 

Rye/ air 

351.12 

Red  fair 

236.05 

Sum  of  Squared 

Deviations 

0.40 

Table  4.4  DBM  Three  Layered  Subcompartment  Vapor  Model  Estimated 
Parameters 


DBM  Three  Layered  Subcompartment  Vapor  Model 


Figure  4.3  DBM  Three  Layered  Subcompartment  Vapor  Model  Results.  This 
graph  shows  how  the  model’s  predictions  compare  to  the  experimen¬ 
tal  data. 
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4.24  Parallel  Subcompartment  Vapor  Model.  The  unknown  subcompart¬ 
ment  permeability  constants  and  partition  coefficients  of  the  parallel  subcompart¬ 
ment  model  were  estimated  and  the  optimized  parameter  values  can  be  seen  in 
Table  4.5.  The  parallel  sub  compartment  model  made  better  predictions  than  the 
homogeneous  model  and  the  layered  sub  compartment  models.  The  parallel  sub¬ 
compartment  model  resulted  in  a  sum  of  squared  deviations  of  0.37  which  is  a  33 
percent  improvement  over  the  homogeneous  model.  Figure  4.4  shows  how  well  the 
predicted  blood  concentration  values  of  the  parallel  subcompartment  model  match 
the  experimental  data. 


Unknown 

Optimized 

Parameters 

Values 

Pfo 

94.04 

Pcd 

0.38 

Rfo/ air 

325.08 

Red  fair 

265.59 

Sum  of  Squared 

Deviations 

0.37 

Table  4.5  DBM  Parallel  Subcompartment  Vapor  Model  Estimated  Parameters 

4.2.5  Parallel-Layered  Subcompartment  Vapor  Model.  The  unknown  sub¬ 
compartment  permeability  constants  and  partition  coefficients  of  the  parallel-layered 
subcompartment  model  were  estimated  and  the  optimized  parameter  values  can  be 
seen  in  Table  4.6.  The  parallel-layered  sub  compartment  model  made  better  blood 
concentration  predictions  than  the  homogeneous  model,  the  layered  subcompartment 
models,  and  the  parallel  subcompartment  model.  The  parallel-layered  subcompart¬ 
ment  model  resulted  in  a  sum  of  squared  deviations  of  0.30  which  is  a  45  percent 
improvement  over  the  homogeneous  model.  Figure  4.5  shows  how  well  the  predicted 
concentration  values  of  the  parallel-layered  subcompartment  model  match  the  ex¬ 
perimental  data. 
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DBM  Parallel  Subcompartment  Vapor  Model 


Figure  4.4  DBM  Parallel  Subcompartment  Vapor  Model  Results.  This  graph  shows 
how  the  model’s  predictions  compare  to  the  experimental  data. 


Unknown 

Optimized 

Parameters 

Values 

Pfo 

94.04 

Psc 

0.46 

Pcd 

2.28 

of  air 

325.08 

R  sc/ air 

649.15 

Red  j  air 

257.93 

Sum  of  Squared 

Deviations 

0.30 

Table  4.6  DBM  Parallel-Layered  Subcompartment  Vapor  Model  Estimated 
Parameters 
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Figure  4.5  DBM  Parallel-Layered  Sub  compartment  Vapor  Model  Results.  This 
graph  shows  how  the  model’s  predictions  compare  to  the  experimental 
data. 
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4.3  Sensitivity  Analysis  for  Vapor  Models 

We  conducted  sensitivity  analysis  on  the  parameters  of  the  skin  compartments 
of  each  model  to  determine  the  impact  they  had  on  the  predicted  blood  concentration 
(Tables  4.7-4.11).  We  looked  at  the  percent  change  in  predicted  blood  concentra¬ 
tions  at  one-half  hour  and  at  four  hours  for  each  change  in  the  parameters  for  the 
500  ppm  and  10,000  ppm  exposures.  In  each  of  the  models,  the  predicted  blood 
concentration  was  the  most  sensitive  to  changes  in  one  of  the  subcompartment  per- 
meabihty  constants.  The  blood  concentration  predictions  were  much  less  sensitive 
to  the  parameters  of  the  mass-balance  equations  for  body  compartments  other  than 
the  skin.  The  most  sensitive  parameters  in  the  other  body  compartments  were  blood 
flow  to  the  liver  (Qi),  fat  blood  partition  coefficient  {Rf/b)i  volume  of  the  fat  com¬ 
partment  (V/),  and  blood  flow  to  the  fat  compartment  (Qf),  in  that  order  (data  not 
shown).  Increasing  the  parameters  by  five  percent  had  relatively  the  same  effect  as 
decreasing  the  parameter  by  five  percent,  although  the  direction  was  different.  The 
rest  of  this  discussion  will  pertain  to  the  effect  of  increasing  the  parameters  by  five 
percent.  Most  of  the  parameters  had  a  greater  sensitivity  at  the  higher  concentra¬ 
tion  than  they  did  at  the  lower  concentration  because  metabohsm  is  saturated  at 
the  higher  concentration  and  the  increase  of  chemical  into  the  body  is  more  directly 
reflected  in  predicted  blood  concentration. 

4-3.1  Homogeneous  Vapor  Model.  In  the  homogeneous  model,  the  order 
of  sensitivity  was  permeabihty  (P),  blood  air  partition  coefficient  {Rb/air)j  skin  air 
partition  coefficient  (Rsk/air),  and  flow  to  the  skin  (Qsk)  (Table  4.7).  Increasing  the 
most  sensitive  parameter,  permeabihty,  increased  the  predicted  blood  concentration 
by  nearly  equivalent  percentages  during  a  one-half  hour  exposure  and  during  a  four- 
hour  exposure.  Changes  in  the  blood  air  and  skin  air  partition  coefficients  and  the 
blood  flow  to  the  skin  had  much  more  effect  on  predicted  blood  concentrations  after 
one-half  hour  than  after  four  hours,  because  early  in  the  exposure  the  concentration 
difference  between  skin  and  the  environment  is  much  greater  than  after  four  hours. 
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500 

ppm 

10,000  ppm 

.5  hours 

4  hours 

.5  hours 

4  hours 

Parameter 

+5% 

+5% 

+5% 

+5% 

P 

5.0 

5.0 

6.4 

5.8 

Rh/aiv 

3.3 

0.5 

4.6 

2.6 

Rsk/air 

-3.1 

-0.1 

-4.1 

-1.0 

Q  sk 

2.7 

0.2 

3.5 

0.5 

Table  4.7  Sensitivity  Analysis  for  the  DBM  Homogenous  Vapor  Model.  These  num¬ 
bers  represent  the  percent  change  in  predicted  blood  concentrations  after 
the  parameters  were  varied. 

4.3.2  Two  Layered  Subcompartment  Vapor  Model.  In  the  two  layered 
subcompartment  model,  changing  the  permeability  for  stratum  corneum  (P*c)  had 
the  largest  impact  on  the  predicted  blood  concentration  while  the  composite  dermal 
permeabihty  (Pcd)  had  practically  no  effect  (Table  4.8).  This  is  because  Pcd  is  so 
large  compared  to  Psc  that  a  five  percent  change  in  Pcd  resulted  in  a  very  small 
change  in  the  overall  permeability  of  the  skin  from  Equation  (3.4),  but  a  five  percent 
change  in  Pgc  resulted  in  about  a  five  percent  change  in  the  overall  permeabihty. 
In  contrast,  the  predicted  blood  concentration  was  much  more  sensitive  to  changes 
in  the  composite  dermal  partition  coefficient  (Pcd/air)  than  it  was  to  changes  in 
the  stratum  corneum  partition  coefficient  {Rsc/air)-  If  is  interesting  to  note  that 
in  this  model  the  volume  of  the  composite  dermal  subcompartment  is  about  fifty 
times  that  of  the  stratum  corneum  subcompartment  and  the  value  for  the  stratum 
corneum  air  partition  coefficient  is  almost  one  order  of  magnitude  greater  than  the 
optimized  value  of  the  composite  dermal  air  partition  coefficient.  In  this  model,  the 
partition  coefficients  and  the  flow  to  the  skin  (Qak)  had  more  impact  on  the  blood 
concentration  after  one-half  hour  than  after  four  hours,  because  the  concentration 
difference  between  air  and  skin  was  diminished  after  four  hours. 

4.3.3  Three  Layered  Subcompartment  Vapor  Model.  In  the  three  layered 
subcompartment  model,  like  the  two  layered  subcompartment  model,  changing  the 
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500 

ppm 

10,000  ppm 

.5  hours 

4  hours 

.5  hours 

4  hours 

Parameter 

+5% 

+5% 

+5% 

+5% 

P 

5.0 

5.0 

6.3 

5.7 

Rb/air 

3.3 

0.6 

4.6 

2.6 

Rcd/air 

-2.7 

-0.1 

-3.6 

-0.9 

Q  sk 

2.7 

0.2 

3.4 

0.5 

R  sc /air 

-0.4 

0.0 

-0.5 

-0.1 

Pcd 

0.0 

0.0 

0.0 

0.0 

Table  4.8  Sensitivity  Analysis  for  the  DBM  Two  Layerd  Subcompartment  Vapor 
Model.  These  numbers  represent  the  percent  change  in  predicted  blood 
concentrations  after  the  parameters  were  varied. 

permeabihty  for  stratum  corneum  {Psc)  had  a  much  larger  impact  on  the  predicted 
blood  concentration  than  the  permeabihty  of  the  viable  epidermis  (Pve)  and  the 
composite  dermal  (Pcd)  (Table  4.9).  When  these  parameters  were  estimated,  the 
composite  dermal  and  the  viable  epidermis  permeabihties  were  significantly  larger 
than  the  stratum  corneum  permeabihty.  This  caused  a  five  percent  change  in  Pcd  and 
P„e  to  change  the  overaU  permeabihty  of  the  skin  by  only  a  smaU  percent  through 
Equation  (3.8)  which  resulted  in  a  low  sensitivity  for  Pcd  and  P„e-  Skin  subcom¬ 
partment  air  partition  coefficients  had  an  order  of  sensitivity  of  composite  dermal 
{Rcd/air),  stratum  corneum  (i2*c/air)j  and  then  viable  epidermis  {Rve/air)-  The  com¬ 
posite  dermal  was  quite  a  bit  more  sensitive.  It  appears  that  the  partition  coefficient 
for  the  composite  dermal  subcompartment  had  the  most  effect  on  chemical  transport 
through  the  rest  of  the  skin. 

^.3.4  Parallel  Subcompartment  Vapor  Model.  In  the  parallel  subcompart¬ 
ment  model,  changing  the  permeabihty  of  the  foUicle  subcompartment  (P/o)  had 
the  largest  effect  on  the  predicted  blood  concentrations  (Table  4.10).  The  blood  air 
partition  coefficient  (Rb/air)  had  the  next  largest  impact  which  is  consistent  with  the 
layered  subcompartment  models.  The  permeabihty  of  the  composite  dermal  subcom- 
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500 

ppm 

10,000  ppm 

.5  hours 

4  hours 

.5  hours 

4  hours 

Parameter 

+5% 

+5% 

+5% 

+5% 

P,C 

5.0 

5.0 

6.3 

5.7 

Rb/air 

3.4 

0.6 

4.6 

2.6 

Rcd/air 

-2.6 

-0.1 

-3.4 

-0.9 

Q  sk 

2.7 

0.2 

3.4 

0.5 

R sc /air 

-0.3 

-0.0 

-0.5 

-0.1 

Rye /air 

-0.2 

0.0 

-0.2 

-0.1 

Pve 

0.0 

0.0 

0.0 

0.0 

Ped 

0.0 

0.0 

0.0 

0.0 

Table  4.9  Sensitivity  Analysis  for  the  DBM  Three  Layered  Sub  compartment  Vapor 
Model.  These  numbers  represent  the  percent  change  in  predicted  blood 
concentrations  after  the  parameters  were  varied. 

partment  (Pcd)  had  the  next  largest  impact  which  was  probably  because  it  contains 
the  stratum  corneum.  The  composite  dermal  air  partition  coefficient  {Red/ air)  was 
more  sensitive  than  the  folhcle  air  partition  coefficient  {Rfo/air)  which  was  probably 
due  to  the  small  fractional  area  of  the  follicle  subcompartment  (10“^).  The  flow 
to  the  folhcle  {Qfo)  and  composite  dermal  {Qcd)  subcompartments  was  much  less 
sensitive  than  the  flow  to  the  skin  (Qsfc)  in  the  layered  subcompartment  models. 

J^.S.S  Parallel-Layered  Subcompartment  Vapor  Model.  In  the  parallel- 

layered  subcompartment  model,  hke  in  the  parallel  subcompartment  model,  changing 
the  permeabihty  of  the  folhcle  sub  compartment  (Pfo)  had  the  largest  effect  on  the 
predicted  blood  concentrations  (Table  4.11).  After  the  blood  air  partition  coeffi¬ 
cient  {Rb/air),  the  permeabihty  of  the  stratum  corneum  sub  compartment  (P,c)  had 
the  next  largest  effect  on  the  blood  concentration  predictions.  Since  the  permeabil¬ 
ity  of  the  composite  dermal  subcompartment  {Pad)  is  relatively  large  compared  to 
Pgc,  changing  Ped  by  five  percent  actually  changed  the  overall  permeabihty  of  the 
non-fofficle  skin  subcompartment  by  a  very  small  percent  because  of  Equation  (3.21). 
The  stratum  corneum  air  partition  coefficient  {Rsc/air)  had  relatively  no  impact  when 
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500 

ppm 

10,000 

ppm 

.5  hours 

4  hours 

.5  hours 

4  hours 

Parameter 

+5% 

+5% 

-f5% 

+5% 

Pfo 

3.9 

3.4 

4.9 

4.0 

^hjair 

1.3 

0.7 

1.8 

2.6 

Pcd 

0.9 

1.5 

0.9 

1.7 

Red /air 

-0.7 

-0.1 

-0.9 

-0.8 

Qio 

0.3 

0.2 

0.4 

0.3 

Qcd 

0.4 

0.1 

0.2 

0.2 

Rfo/air 

0.0 

0.0 

-0.1 

0.0 

Table  4.10  Sensitivity  Analysis  for  the  DBM  Parallel  Sub  compartment  Vapor 
Model.  These  numbers  represent  the  percent  change  in  predicted  blood 
concentrations  after  the  parameters  were  varied. 

varied  which  is  consistent  with  the  layered  subcompartment  models.  The  remaining 
parameters  had  the  same  order  of  sensitivity  and  relatively  the  same  magnitudes  of 
sensitivity  as  in  the  parallel  subcompartment  model. 


^-4  Liquid  Model  Predictions 

The  blood  concentration  predictions  from  the  five  hquid  models  were  compared 
to  previously  collected  data  in  which  rats  were  exposed  to  neat  dibromomethane 
(DBM)  hquid  in  a  cell  with  an  area  of  3.14  cm  at  a  dose  of  7434  ppm  for  a  twenty  four 
hour  period  [16].  Although  the  hquid  and  vapor  models  both  use  dibromomethane 
(DBM),  they  should  not  be  compared  because  of  the  differences  in  some  of  the 
parameter  values  for  DBM.  These  parameter  differences  are  a  result  of  different 
researchers  measuring  the  parameters  and  a  ten  year  gap  in  the  development  of  the 
models. 

The  subcompartment  models  were  optimized  using  the  technique  described  in 
Section  3.3,  however  we  weighted  the  last  data  point  very  low  (100  times  lower) 
compared  to  the  first  seven  data  points.  We  did  this  so  the  optimization  technique 
would  focus  on  the  early  data  points  since  those  are  the  times  we  are  most  interested 
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500 

ppm 

10,000  ppm 

.5  hours 

4  hours 

.5  hours 

4  hours 

Parameter 

+5% 

+5% 

+5% 

+5% 

Pfo 

3.8 

3.2 

4.7 

3.7 

^bjair 

1.3 

0.7 

1.9 

2.6 

P.C 

1.0 

1.7 

1.1 

1.9 

^cd/ air 

-0.8 

-0.1 

-1.0 

-0.8 

Qfo 

0.3 

0.2 

0.4 

0.2 

Qcd 

0.5 

0.1 

0.3 

0.2 

Rfo/air 

0.2 

0.2 

0.3 

0.3 

Rscjair 

0.0 

0.0 

-0.1 

0.0 

Pcd 

0.0 

0.0 

0.0 

0.0 

Table  4.11  Sensitivity  Analysis  for  the  DBM  Parallel-Layered  Subcompartment  Va¬ 
por  Model.  These  numbers  represent  the  percent  change  in  predicted 
blood  concentrations  after  the  parameters  were  varied. 

in.  Also,  we  are  not  sure  if  the  chemical  physically  changes  the  skin  after  a  long 
exposure.  Figures  4.6-4.10  show  how  each  model  compares  to  the  experimental  data 
individually  and  Figure  4.11  shows  all  five  of  the  hquid  models  on  the  same  graph 
to  illustrate  the  relative  improvements. 

4-4  1  Homogeneous  Liquid  Model.  The  homogeneous  hquid  model  skin 
compartment  is  set  up  to  represent  the  stratum  corneum  only,  as  opposed  to  the 
homogeneous  vapor  model  which  has  a  skin  compartment  representing  the  whole 
skin.  The  parameters  in  the  homogeneous  hquid  model  present  a  relatively  thin  skin 
compartment  with  the  properties  of  the  upper  layer  of  the  skin. 

The  parameters  for  the  homogeneous  model  were  determined  experimentaUy 
or  found  in  the  hterature  (Tables  4.12  and  4.13)  [16].  It  should  be  noted  that  the 
permeabihty  of  the  skin  in  this  model  was  optimized  using  an  equal  weight  on  aU 
of  the  data  points  (Pgk  =  .003).  We  also  optimized  the  permeabihty  using  the 
same  technique  that  was  used  for  optimizing  parameters  in  the  subcompartment 
models  (lower  weight  on  the  last  data  point)  and  obtained  a  different  permeabihty 
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(P,fc  =  .0032).  Figure  4.6  shows  how  the  predictions  of  the  previously  developed 
homogeneous  model  compare  to  the  experimental  data  using  both  values  for  the 
permeabihty.  For  the  two  layered  subcompartment  model,  we  used  the  original 
permeabiUty  {Pgk  =  .003)  as  the  permeability  of  the  stratum  corneum  and  after 
optimizing  the  permeability  of  the  composite  dermal  subcompartment,  we  calculated 
the  overall  skin  permeability  using  Equation  (3.4)  which  we  then  used  for  the  other 
subcompartment  models.  For  this  reason,  we  wiU  compare  the  subcompartment 
models  to  the  original  homogeneous  model  which  uses  Psk  =  .003.  The  original 
homogeneous  model  predicted  the  concentrations  in  the  blood  reasonably  well  having 
a  sum  of  squared  deviations  of  0.098  when  compared  to  the  experimental  data. 


Compartment 

Tissue/air  partition 
coefficient 

Blood  flow 
(  %  cardiac  output) 

Volume 

(  %  body  weight) 

Skin 

120.00 

5 

10 

Rapidly  Perfused 

0.9 

56 

5 

Slowly  Perfused 

0.546 

10 

65 

Fat 

10.8 

9 

7 

Liver 

0.918 

20 

4 

Blood 

74.1 

N/A 

N/A 

Table  4.12  DBM  Parameters  for  the  tissue  compartments  of  the  homogenous  liquid 
model 


Constants 

Values 

Vmax  (mg/hr/kg) 

12.5 

Km  (mg/hter) 

.36 

K,  (hr-' -kg-') 

.557 

Kp  (cm /hr) 

0.003 

Table  4.13  DBM  Constants  for  the  homogenous  Hquid  model 


>^.4-2  Two  Layered  Subcompartment  Liquid  Model.  The  unknown  sub¬ 
compartment  permeabihty  constants  and  partition  coefficients  of  the  two  layered 
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DBM  Homogenous  Liquid  Model 


time(hours) 

Figure  4.6  DBM  Homogeneous  Liquid  Model  Results.  This  graph  shows  how  the 
model’s  predictions  compare  to  the  experimental  data. 

subcompartment  model  were  estimated  and  the  optimized  parameter  values  can  be 
seen  in  Table  4.14.  The  two  layered  sub  compartment  model  matched  the  data  better 
than  the  homogeneous  model,  yielding  a  sum  of  squared  deviations  of  0.063  when 
compared  to  the  experimental  data.  The  two  layered  subcompartment  model’s  sum 
of  squared  deviations  is  a  36  percent  improvement  in  accuracy  over  the  homogeneous 
model.  Figure  4.7  shows  how  well  the  predicted  values  of  the  two  layered  subcom¬ 
partment  model  match  the  experimental  data.  The  two  layered  subcompartment 
model  predicts  the  concentration  of  the  chemical  in  the  blood  more  accurately  than 
the  homogeneous  model,  because  it  is  a  more  physiologically  detailed  representation 
of  the  skin  and  it  has  one  more  parameters  that  was  optimized. 

Three  Layered  Subcompartment  Liquid  Model  The  unknown  sub¬ 
compartment  permeability  constants  and  partition  coefficients  of  the  three  layered 
sub  compartment  model  were  estimated  and  the  optimized  parameter  values  can 


Unknown 

Parameters 

Optimized 

Values 

Psc 

0.003 

Pcd 

5.00 

^ SC /air 

120.00 

Red/ air 

66.96 

Sum  of  Squared 
Deviations 

0.063 

Table  4.14  DBM  Two  Layered  Subcompartment  Liquid  Model  Estimated 
Parameters 


DBM  Two  Layered  Subcompartment  Liquid  Model 


Figure  4.7  DBM  Two  Layered  Subcompartment  Liquid  Model  Results.  This  graph 
shows  how  the  model’s  predictions  compare  to  the  experimental  data. 
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be  seen  in  Table  4.15.  The  three  layered  subcompartment  model  yielded  about  the 
same  results  as  the  two  layered  subcompartment  model  resulting  in  a  sum  of  squared 
deviations  of  0.060  when  compared  to  the  experimental  data.  The  three  layered  sub¬ 
compartment  model’s  sum  of  squared  deviations  is  a  39  percent  improvement  in 
accuracy  over  the  homogeneous  model.  Figure  4.8  shows  how  well  the  predicted 
concentration  values  of  the  three  layered  subcompartment  model  match  the  experi¬ 
mental  data.  The  three  layered  subcompartment  model  yields  more  accurate  blood 
concentration  predictions,  as  well  as,  a  more  detailed  representation  of  the  structure 
of  the  skin. 


Unknown 

Optimized 

Parameters 

Values 

P.C 

0.003 

Pve 

7.50 

Pcd 

14.99 

^  SC  fair 

120.00 

^ve/ air 

102.00 

Red/ air 

65.50 

Sum  of  Squared 

Deviations 

0.060 

Table  4.15  DBM  Three  Layered  Subcompartment  Liquid  Model  Estimated 
Parameters 


Parallel  Subcompartment  Liquid  Model.  The  unknown  subcompart¬ 
ment  permeabihty  constants  and  partition  coefficients  of  the  parallel  subcompart¬ 
ment  model  were  optimized  and  the  optimized  parameter  values  can  be  seen  in 
Table  4.16.  The  parallel  subcompartment  model  made  better  predictions  than  the 
homogeneous  model  and  the  layered  subcompartment  models.  The  parallel  sub¬ 
compartment  model  resulted  in  a  sum  of  squared  deviations  of  0.054  which  is  a  45 
percent  improvement  over  the  homogeneous  model.  Figure  4.9  shows  how  well  the 
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time  (hours) 


Figure  4.8  DBM  Three  Layered  Subcompartment  Liquid  Model  Results.  This 
graph  shows  how  the  model’s  predictions  compare  to  the  experimen¬ 
tal  data. 


predicted  blood  concentration  values  of  the  parallel  subcompartment  model  match 
the  experimental  data. 


Unknown 

Optimized 

Parameters 

Values 

Pfo 

0.07 

Pcd 

0.002 

Rfo/air 

271.67 

Red  1  air 

66.58 

Sum  of  Squared 

Deviations 

0.054 

Table  4.16  DBM  Parallel  Subcompartment  Liquid  Model  Estimated  Parameters 


DBM  Parallel  Subcompartment  Liquid  Model 


Figure  4.9  DBM  Parallel  Subcompartment  Liquid  Model  Results.  This  graph 
shows  how  the  model’s  predictions  compare  to  the  experimental  data. 

>{.4 -5  Parallel-Layered  Subcompartment  Liquid  Model.  The  unknown  sub¬ 
compartment  permeability  constants  and  partition  coefficients  of  the  parallel-layered 
subcompartment  model  were  optimized  and  the  optimized  parameter  values  can  be 
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seen  in  Table  4.17.  The  parallel-layered  sub  compartment  model  made  better  blood 
concentration  predictions  than  the  homogeneous  model  and  the  layered  subcom¬ 
partment  models,  but  not  quite  as  well  as  the  parallel  sub  compartment  model.  The 
parallel-layered  subcompartment  model  resulted  in  a  sum  of  squared  deviations  of 
0.058  which  is  a  41  percent  improvement  over  the  homogeneous  model.  Figure  4.10 
shows  how  well  the  predicted  blood  concentration  values  of  the  parallel-layered  sub¬ 
compartment  model  match  the  experimental  data. 


Unknown 

Optimized 

Parameters 

Values 

Pfo 

0.07 

Psc 

0.002 

Pcd 

5.13 

Rfo/air 

271.67 

Rsc/aiv 

146.47 

^cdfair 

64.98 

Sum  of  Squared 

Deviations 

0.058 

Table  4.17  DBM  Parallel-Layered  Subcompartment  Liquid  Model  Estimated 
Parameters 


^.5  Sensitivity  Analysis  for  Liquid  Models 

We  conducted  sensitivity  analysis  on  the  parameters  of  the  skin  compartments 
of  each  model  to  determine  the  impact  they  have  on  the  predicted  blood  concen¬ 
tration  (Tables  4.18-4.22).  We  looked  at  the  percent  change  in  predicted  blood 
concentrations  at  one-half  hour  and  at  twenty  four  hours  for  each  change  in  the 
parameters.  In  each  of  the  models,  the  predicted  blood  concentration  was  the  most 
sensitive  to  changes  in  one  of  the  subcompartment  permeabihty  constants.  The  pa¬ 
rameters  of  the  mass-balance  equations  for  body  compartments  other  than  the  skin 
had  relatively  the  same  sensitivity  as  they  did  in  the  vapor  model.  Increasing  the 
paxameters  by  five  percent  had  relatively  the  same  effect  as  decreasing  the  parameter 
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DBM  Comparison  of  All  Five  Liquid  Models 


Figure  4.11  DBM  Comparison  of  the  Results  of  All  Five  Liquid  Models.  This  graph 
shows  how  the  model’s  predictions  compare  to  the  experimental  data. 

by  live  percent,  although  the  direction  was  different.  The  rest  of  this  discussion  will 
pertain  to  the  effect  of  increasing  the  parameters  by  five  percent. 

4.5.1  Homogeneous  Liquid  Model.  In  the  homogeneous  model,  the  order 
of  sensitivity  was  permeabihty  (P),  blood  air  partition  coefficient  {Rb/air)j  skin  air 
partition  coefficient  (Rsk/air),  and  flow  to  the  skin  (Q^k)  (Table  4.18).  Increasing 
the  permeabihty,  the  skin  air  pajtition  coefficient,  and  the  flow  to  the  skin  increased 
the  predicted  blood  concentration  by  nearly  equivalent  percentages  during  a  one-half 
hour  exposure  and  during  a  twenty  four  hour  exposure.  Opposite  of  the  vapor  model, 
changes  in  the  blood  air  partition  coefficient  had  a  greater  effect  after  twenty  four 
hours  than  after  one-half  hour. 

4.5.2  Two  Layered  Suhcompartment  Liquid  Model.  In  the  two  layered 
subcompartment  model,  changing  the  permeabihty  for  stratum  corneum  (Psc)  had 
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Parameter 

.5  hours 
-f5% 

24  hours 
-h5% 

p 

5.6 

5.6 

Rb/air 

1.3 

5.1 

Rsk/air 

-0.8 

-0.8 

Q  sk 

0.8 

0.8 

Table  4.18  Sensitivity  Analysis  for  the  DBM  Homogenous  Liquid  Model.  These 
numbers  represent  the  percent  change  in  predicted  blood  concentrations 
after  the  parameters  were  varied. 

the  largest  impact  on  the  predicted  blood  concentration  while  the  composite  dermal 
permeability  (Pcd)  had  practically  no  effect  (Table  4.19).  This  was  because  Pcd  is  so 
large  compared  to  Pgc  that  a  five  percent  change  in  Pcd  resulted  in  a  very  small  change 
in  the  overall  permeability  of  the  skin  from  Equation  (3.4),  but  a  five  percent  change 
in  Psc  resulted  in  about  a  five  percent  change  in  the  overall  permeability.  In  contrast, 
the  predicted  blood  concentration  was  more  sensitive  to  changes  in  the  composite 
dermal  partition  coefficient  (Rcd/air)  than  it  was  to  changes  in  the  stratum  corneum 
partition  coefficient  (Rsc/air)-  As  in  the  homogeneous  model,  the  blood  air  partition 
coefficient  (Rb/air)  had  a  greater  effect  after  twenty  four  hours  than  after  one-half 
hour.  In  this  model,  the  composite  dermal  partition  coefficient  and  the  blood  flow 
to  the  skin  (Qnk)  had  more  impact  on  the  blood  concentration  after  one-half  hour 
than  after  twenty  four  hours,  because  the  concentration  difference  between  air  and 
skin  was  diminished  after  twenty  four  hours. 

4-5.3  Three  Layered  Subcompartment  Liquid  Model.  In  the  three  layered 
subcompartment  model,  like  the  two  layered  subcompartment  model,  changing  the 
permeability  for  stratum  corneum  (Psc)  had  a  much  larger  impact  on  the  predicted 
blood  concentration  than  the  permeability  of  the  viable  epidermis  (P^e)  and  the 
composite  dermal  (Pcd)  (Table  4.20).  When  these  parameters  were  estimated,  the 
composite  dermal  and  the  viable  epidermis  permeabilities  were  significantly  larger 
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.5  hours 

24  hours 

Parameter 

+5% 

-f5% 

P.C 

5.8 

5.8 

Rb/air 

1.4 

4.8 

Q  sk 

0.9 

0.5 

Rcd/air 

-0.5 

0.1 

Pcd 

0.1 

0.1 

Rsc/air 

0.0 

0.0 

Table  4.19  Sensitivity  Analysis  for  the  DBM  Two  Layerd  Subcompartment  Liquid 
Model.  These  numbers  represent  the  percent  change  in  predicted  blood 
concentrations  after  the  parameters  were  varied. 

than  the  stratum  corneum  permeabihty.  This  caused  a  five  percent  change  in  Pcd  and 
P„e  to  change  the  overall  permeabihty  of  the  skin  by  only  a  small  percent  through 
Equation  (3.8)  which  resulted  in  a  very  low  sensitivity  for  Pcd  and  P„e.  On  the  other 
hand,  the  composite  dermal  air  partition  coefficient  {Rcd/air)  was  the  only  subcom¬ 
partment  partition  coefficient  to  have  any  effect  on  the  predicted  blood  concentration 
when  changed.  The  effect  of  changing  the  blood  flow  to  the  skin  {Qsk)  was  greater 
after  one-half  hour  than  after  twenty  four  hours.  As  in  all  of  the  hquid  models,  the 
blood  air  partition  coefficient  had  a  larger  effect  after  one-half  hour  than  after  twenty 
four  hours. 

4-5.4  Parallel  Subcompartment  Liquid  Model.  In  the  parallel  subcompart¬ 
ment  model,  changing  the  permeability  of  the  composite  dermal  subcompartment 
(Pcd)  had  the  largest  effect  on  the  predicted  blood  concentrations  probably  because 
it  contains  the  stratum  corneum  (Table  4.10).  The  blood  air  partition  coefficient 
(Rb/air)  had  the  next  largest  impact  which  is  consistent  with  the  layered  subcom¬ 
partment  models.  The  permeabihty  of  the  folhcle  subcompartment  (P/o)  had  the 
next  largest  impact  which  is  probably  because  it  acts  as  its  own  small  skin  com¬ 
partment  having  the  chemical  enter  from  the  surface  and  leave  to  the  blood.  The 
composite  dermal  air  partition  coefficient  (Rcd/air)  was  more  sensitive  than  the  foUi- 
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.5  hours 

24  hours 

Parameter 

-F5% 

-f-5% 

Psc 

5.8 

5.8 

Rb/ air 

1.4 

4.8 

Q  sk 

0.9 

0.4 

^cdfaiT 

-0.5 

0.0 

P 

ve 

0.0 

0.0 

Rye /air 

0.0 

0.0 

Rsc/air 

0.0 

0.0 

Pcd 

0.0 

0.0 

Table  4.20  Sensitivity  Analysis  for  the  DBM  Three  Layered  Subcompartment  Liq¬ 
uid  Model.  These  numbers  represent  the  percent  change  in  predicted 
blood  concentrations  after  the  parameters  were  varied. 

cle  air  partition  coefficient  (ii/o/air)  which  was  probably  due  to  the  small  fractional 
area  of  the  foUicle  subcompartment  (10“^).  The  flow  to  the  composite  dermal  {Qcd) 
subcompartment  had  about  the  same  effect  as  the  flow  to  the  composite  dermal 
subcompartments  of  the  layered  models,  but  the  flow  to  the  folhcle  {Qfo)  subcom¬ 
partment  had  a  somewhat  lower  sensitivity. 


.5  hours 

24  hours 

Parameter 

4-5% 

-b5% 

Pcd 

4.4 

4.5 

Rh/air 

1.4 

4.8 

Pfo 

1.5 

1.4 

Qcd 

0.9 

0.3 

Red /air 

-0.6 

0.0 

Qfo 

0.1 

0.1 

Rfo/air 

-0.1 

-0.1 

Table  4.21  Sensitivity  Analysis  for  the  DBM  Parallel  Subcompartment  Liquid 
Model.  These  numbers  represent  the  percent  change  in  predicted  blood 
concentrations  after  the  parameters  were  varied. 
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^.5.5  Parallel-Layered  Subcompartment  Liquid  Model.  In  the  parallel- 
layered  subcompartment  model,  changing  the  permeabihty  of  the  stratum  corneum 
sub  compartment  (P«c)  had  the  largest  effect  on  the  predicted  blood  concentrations 
(Table  4.22).  After  the  blood  air  partition  coefficient  (Rb/air),  the  permeability  of  the 
folhcle  subcompartment  (P/o)  had  the  next  largest  effect  on  the  blood  concentration 
predictions.  Since  the  permeabihty  of  the  composite  dermal  subcompartment  (Pcd) 
is  relatively  large  compared  to  P^c,  changing  Pcd  by  five  percent  actually  changed 
the  overall  permeabihty  of  the  non-foUicle  skin  subcompartment  by  a  very  smaU  per¬ 
cent  because  of  Equation  (3.21).  The  stratum  corneum  air  {Rsc/air)  and  folhcle  air 
(Rfo/air)  partition  coefficients  had  relatively  no  impact  when  varied.  The  remaining 
parameters  had  the  same  order  of  sensitivity  and  relatively  the  same  magnitudes  of 
sensitivity  as  in  the  paraUel  subcompartment  model. 


.5  hours 

24  hours 

Parameter 

+5% 

4-5% 

Psc 

4.3 

4.4 

Rb/air 

1.4 

4.8 

Pfo 

1.5 

1.4 

Qcd 

0.9 

0.3 

Red /air 

-0.5 

0.0 

Qfo 

0.1 

0.1 

Rsc/air 

0.0 

0.0 

Rfo/air 

0.0 

0.0 

Pcd 

0.0 

0.0 

Table  4.22  Sensitivity  Analysis  for  the  DBM  ParaUel-Layered  Subcompartment 
Liquid  Model.  These  numbers  represent  the  percent  change  in  predicted 
blood  concentrations  after  the  parameters  were  varied. 


4.6  Summary 

The  vapor  and  hquid  models  were  used  with  parameters  from  studies  involv¬ 
ing  dibromomethane  (DBM).  The  results  from  all  of  the  models  were  compared,  by 
graphical  plots  and  calculating  the  sura  of  squared  deviations,  to  the  experimental 
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data  that  was  previously  collected  in  the  laboratory  .  The  sub  compartment  models 
were  an  improvement  over  the  homogeneous  models.  We  cannot  make  any  compar¬ 
isons  between  the  vapor  and  hquid  models  in  this  research  because  different  values 
for  the  DBM  parameters  were  used  when  the  two  original  models  were  developed. 
Sensitivity  analysis  shows  how  changing  the  skin  subcompartment  parameters  effect 
the  blood  concentration  predictions. 
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V.  Summary  and  Conclusions 


5. 1  Summary 

The  goal  of  this  research  was  to  develop  a  dermal  PBPK  model  with  multiple 
dermal  subcompartments  that  improved  prediction  of  mammalian  blood  concentra¬ 
tions  of  a  chemical  following  dermal  vapor  or  neat  hquid  exposure  and  provide  the 
capability  to  predict  the  concentration  of  chemicals  in  distinct  parts  of  the  skin.  This 
research  resulted  in  the  development  of  four  new  models  which  were  adjusted  to  han¬ 
dle  vapor  and  neat  liquid  exposures.  A  two  layered  subcompartment,  three  layered 
subcompartment,  parallel  subcompartment,  and  parallel-layered  subcompartment 
model  was  developed.  These  new  models  were  derived  from  previously  developed 
homogeneous  models.  The  predictions  of  aU  of  the  models  were  compared  to  the 
experimental  data  in  which  rats  were  exposed  to  dibromomethane  vapor  and  neat 
Hquid  dibromomethane  in  doses  and  time  durations  matching  those  used  in  the  mod¬ 
els.  Sensitivity  analysis  of  the  models’  parameters  showed  which  parameters  are  the 
most  important  to  identify. 

5. 2  Conclusions 

Each  new  model  in  turn  predicted  blood  concentrations  of  dibromomethane 
more  accurately  than  the  original  homogeneous  models  especially  at  the  beginning 
of  the  exposures.  The  early  time  periods,  before  the  chemical  concentration  in  the 
blood  reaches  steady-state,  are  the  most  important  to  model  because  most  realistic 
exposures  will  not  last  for  extended  periods  of  time.  The  new  models  also  represent 
the  structure  of  the  skin  in  more  detail,  making  the  new  models  more  physiologically 
correct.  The  improvement  in  prediction  of  the  three  subcompartment  models  (three 
layered  or  parallel-layered)  over  the  two  subcompartment  models  (two  layered  or  par¬ 
allel)  was  not  nearly  as  great  as  the  improvement  of  the  two  subcompartment  models 
over  the  homogeneous  models.  The  simplest  model  that  works  should  be  used,  and 
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this  argues  for  one  of  the  two  subcompartment  models  (two  layered  or  parallel). 
Choice  of  model  really  depends  on  the  chemical  being  modeled,  the  purpose  of  the 
model,  and  whether  the  additional  model  parameters  for  the  three  sub  compartment 
models  can  readily  be  measured  or  estimated. 

Sensitivity  analysis  provided  insight  into  the  most  important  parameters  to 
measure  or  estimate.  For  these  dermal  exposures,  changing  the  skin  subcompartment 
parameters  had  more  impact  on  the  predicted  blood  concentrations  than  changing 
the  parameters  in  the  other  body  compartments  of  the  model.  In  aU  four  models, 
one  of  the  permeability  constants  was  the  most  sensitive  to  change,  followed  by  the 
blood  air  partition  coefficient.  These  parameters  should  be  concentrated  on  the  most 
when  measuring  or  estimating  parameters  for  these  types  of  models. 

Prediction  of  therapeutic  or  toxic  events  in  the  skin  requires  an  understand- 
ing  of  the  distribution  of  chemicals  in  the  skin.  Expanded  biologically-based  skin 
models  such  as  these  could  be  used  as  the  basis  for  predicting  efficacy  of  a  drug 
which  might  act  differently  in  specific  parts  of  the  skin.  These  models  might  also  be 
used  to  predict  deleterious  dermal  effects  such  as  dermatitis  or  contact  sensitization. 
Mechanistic  understanding  of  the  therapeutic  or  toxic  events,  as  well  as,  their  con¬ 
centration  dependence  will  be  required,  but  physiologically-based  pharmacokinetic 
models  such  as  these  would  be  the  first  step.  Validation  of  model  predictions  for 
chemical  concentration  in  the  skin  subcompartments  is  necessary  and  should  be  ac¬ 
complished  for  several  chemicals  with  different  properties  and  at  different  exposure 
doses  and  durations  before  this  model  can  be  used  mechanistically. 

Permeability  constants  and  partition  coefficients  for  the  new  subcompartments 
were  estimated  or  optimized  in  these  new  models.  It  is  essential  to  be  able  to 
experimentally  determine  these  parameters  if  the  ability  of  a  model  to  extrapolate 
to  humans  is  to  be  maintained.  Measuring  the  partition  coefficients  for  the  skin 
subcompartments  involves  technical  challenges  which  probably  can  be  overcome, 
however  currently  we  can  constrain  the  sub  compartment  partition  coefficients  by 
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the  overall  skin  air  partition  coefficient  using  the  relationship  of  Equation  (2.8). 
Measurement  of  permeabihty  constants  for  each  subcompartment  is  probably  not 
possible,  however  constraining  the  estimates  by  the  overall  permeabihty  constant 
according  to  the  series  relationships  shown  in  Equations  (2.5)  and  (2.6)  will  not  be 
a  problem  since  the  overall  permeability  can  be  estimated  using  the  model. 

Biologically-based  models  have  tremendous  potential  to  provide  species,  dose, 
and  duration  extrapolations  of  laboratory  experiments.  Success  in  these  endeavors 
could  have  a  big  impact  on  pharmacology  and  toxicology.  The  strength  of  these 
models  is  that  they  can  be  predictive,  as  well  as,  descriptive.  But,  the  keys  to 
accurate  prediction  are  a  sound  biological  basis,  and  including  only  parameters  in  the 
model  which  can  be  measured  (or  estimated  if  absolutely  necessary)  in  laboratory 
animals  and  in  humans.  These  models  with  biologically-based  subcompartments 
provide  improved  fit  to  the  data  with  the  potential  for  species  extrapolation  after 
appropriate  measurements  of  human  parameters. 

5.3  Recommendations 

A  vapor  exposure  model  using  the  same  chemical  parameters  presented  in  the 
neat  hquid  homogeneous  model  would  allow  a  comparison  between  the  improvements 
of  the  subcompartment  models  for  different  types  of  exposure.  The  apphcation  of 
the  models  developed  in  this  research  to  chemicals  with  different  properties  could 
give  some  insight  on  how  chemicals  with  certain  properties  penetrate  the  skin. 

The  technique  of  using  subcompartments  to  represent  tissue  compartments, 
which  are  identified  in  a  model,  can  be  used  in  future  research  to  more  accurately  rep¬ 
resent  other  tissue  compartments  and  the  complex  functions  that  take  place  within 
them.  For  example,  the  process  of  metabohsm  and  chemical  binding  could  more 
accurately  be  represented. 
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Appendix  A.  Symbols  for  Equations 

The  following  is  a  complete  list  of  the  symbols  used  in  the  equations  in  this 
thesis. 

D  =  Average  membrane  diffusion  coefficient  (distance^  time~^) 

X  =  Distance 

Js  =  Steady-state  flux  of  solute  (moles  distance”^  time~^) 

S  =  membrane  thickness  (cm) 

AC's  =  Concentration  difference  of  solute  across  membrane  (moles  cm""^) 

Km  =  (Solute  sorbed  per  cc  of  tissue) /(Solute  in  solution  per  cc  of  solvent) 
kp  =  Permeability  constant  (distance/time) 
a  =  arterial 
A  =  area  (distance^) 
al  =  alveolar 
b  =  blood 

C  —  concentration  (mass/volume) 
c  =  cardiac 
/  =  fat 

i  =  tissue  compartment  i 

K  =  Michaehs  constant  -  metabolism  (mass/volume) 
kf  =  first-order  metabohc  rate  constant  (time“^) 

I  =  liver 

P  —  permeability  (distance/time) 
p  =  pulmonary 
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Q  =  flow  (volume/time) 

R  =  partition  coefficient  (ratio  of  concentrations) 
rp  =  rapidly  perfused 
sfc  =  surface 
sk  =  skin 

sp  =  slowly  perfused 
t  =  time 

V  =  volume 

V  =  venous 

Vmax  =  maximum  reaction  velocity  (mass/time) 

sc  =  stratum  corneum 

cd  —  composite  dermal 

ve  =  viable  epidermis 

fo  =  Follicle 

Iq  =  Liquid 

SSD  =  Sum  of  squared  deviations 
n  =  Number  of  experimental  data  points 
TABS  =  Total  amount  absorbed 
AM  =  Amount  (mass) 
ex  =  Exhaled 
met  =  Metabolized 
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